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ABSTRACT 



The Department of Aeronautics and Astronautics at the Naval Postgraduate School (NPS) is 
expanding its helicopter research capabilities in order to facilitate present and future research 
demands. The rapidly changing needs have already out paced available assets. Therefore it was 
necessary to design and develop a new remotely piloted helicopter (RPH) that would meet present 
needs, NOTAR and HHC, and be flexible enough to meet future needs. The research efforts 
encompassed by this thesis are defining the present needs, investigating what type/size of RPH would 
fulfill these needs, procuring this asset, and analyzing its capabilities. Based on a defined payload, 
helicopter trends are analyzed to determine an estimate of the overall RPH size (gross weight) and 
engine size required. A preliminary design process validates these figures. Choosing to procure an 
RPH instead of building one, a detailed performance analysis is conducted on the main rotor system. 
This analysis includes blade vibration analysis, retreating blade stall analysis, and power required 
analysis. Modification of the RPH’s main rotor hub, drive train, and landing gear are studied and 
recommendations presented. This research effort is a continuation of a long- term program to provide 
NPS with robust assets to support present and future rotorcraft research efforts. 
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I. 



INTRODUCTION 



Helicopter research in the Department of Aeronautics and 
Astronautics at the Naval Postgraduate School (NPS) is 
rapidly accelerating. This fast pace has created a need for 
expanding the capabilities of the department's Remotely 
Piloted Vehicles (RPV) . Valuable research has been 
accomplished through the employment of the existing RPVs . 

The GMP Legend, a commercially produced radio controlled 
(RC) helicopter, was used to produce baseline vibration 
analysis and to validate measurement techniques. The 
results from this work will be of great value in subsequent 
Higher Harmonic Control (HHC) research. The Bruiser, a 
limited production, 20 pound payload RPV developed by 
Pacific RPV, was used in a shake test to obtain an airframe 
modal analysis. 

The preceding research has been of great value to the 
department, but in order to advance beyond proof of concept 
and into scale model analysis, a new helicopter was 
required. Incorporated in the requirements were a 
helicopter with a tip speed and chord which produced a 
Reynolds number that was within an acceptable range for 
comparison with a full scale helicopter; a helicopter of a 
size large enough to carry a No-Tail -Rotor (NOTAR) tail boom 
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that would also be testable in the quarter scale (14 by 22 
foot) wind tunnel at NASA Langley; and a helicopter with 
enough flexibility to be adaptable to future research needs. 
This impetus developed the challenge of designing a suitable 
RPV for the department. 

A critical amount of background knowledge and direction 
was obtained on a research trip to the Aerostructures 
Directorate at NASA Langley. They were heavily involved in 
the design, manufacture and test of their own Free-Flight 
Rotorcraft Research Vehicle (FFRRV) , which was directly 
along the same lines of interest of NPS. Their RPV was 
twice the size and considerably more complex than that 
desired by NPS. The head of the directorate, Arthur E. 
Phelps III, was invaluable in passing on their corporate 
knowledge, saving untold hours in achieving the NPS goal. 
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II. 



BACKGROUND 



A. HIGHER HARMONIC CONTROL (HHC) 

The control of vibrations has always been of great 
concern to both the helicopter designer and to the 
helicopter pilot. It has been a continuing source of 
agitation throughout the years and a focus of enormous 
amounts of research assets. The current means by which the 
vibrations are reduced are through passive devices which 
either isolate the source of vibration (isolators) or 
diffuse the vibration level (absorbers) . These vibration 
absorption mechanisms are usually restricted to a narrow 
scope of flight conditions and vibratory frequencies. The 
use of HHC, a relatively new technology, is an active 
vibration reduction device, vice the passive ones just 
mentioned. HHC functions by altering the aerodynamic loads 
on the rotor, and therefore the vibratory forces and moments 
which cause the airframe to vibrate are reduced. [Ref. 1] 

In earlier full scale HHC testing on the Hughes 
Helicopter 0H-6A, it was determined that not only were the 
vibrations successfully reduced, but the performance of the 
helicopter was also improved. Other helicopter companies 
tried to duplicate this resulting improved performance but 
were unsuccessful. Corroborating this performance 
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enhancement locally would have great merit. The RPVs on 
hand did not allow this type of research to be conducted 
because they could not produce Reynolds numbers over the 
main rotor blade that would allow data comparison with the 
data obtained on the 0H-6A. A helicopter with a tip speed 
and chord length that would produce Reynolds numbers on the 
order of two million was required. This shortcoming 
provided some of the drive to obtain a new helicopter RPV. 

B . NOTAR 

Presently, Lt . M. Borno, in conjunction with McDonnell 
Douglas Helicopter Corporation, is conducting thesis 
research which will produce a fully operational scale model 
NOTAR tailboom. The tailboom requirement was that it was to 
be big enough to test in the quarter scale wind tunnel at 
NASA Langley. The current RPVs lacked the size to 
mechanically support and drive a tailboom of this scale. 
This was a second driver behind acquiring a new helicopter 
RPV. 



C. AUTONOMOUS LANDING AND TAKEOFF SYSTEM (ALTOS) 

The United States Navy's growing interest in RPVs and 
their fleet applications has now reached a point which 
requires significant background research. As the Navy leans 
toward fully autonomous RPVs, it lacks experience and 
expertise in the most effective and efficient means to 
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launch and recover these vehicles. Recently, Orion Aviation 
entered into a contract with the United States Navy in 
conjunction with NPS to develop and test five potential 
ALTOS concepts and have them analyzed for merit by students 
at NPS. Once the best ALTOS concept is determined, it will 
be built by Orion Aviation and then will be demonstrated 
using an NPS RPV. The expanded capabilities that the new 
NPS RPV would add to the existing resources would provide 
greater flexibility in the demonstration and validation 
phase of the ALTOS program. 
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III. 



DEVELOPING SCALE MODEL TRENDS 



There is a great deal of documentation for the design of 
full scale helicopters regardless of weight range, but there 
exits little or no documentation for design of a scale model 
helicopter. It was necessary to determine trends for gross 
weight, takeoff weight, disk loading (DL) , rotor radius, 
solidity (cr), and blade loading (BL) from full scale 
helicopters in order to determine what the values should be 
for quarter- scale size. The trends that were developed were 
typically linear with an adjustment of either the x- 
intercept or the y- intercept usually required. The 
adjustment was made assuming the trends held true for all 
weight ranges, but at the lower weights the line had to be 
shifted to accommodate scale sizes. The trends that were 
generated show this assumption to be good. 

The first trend determined was between the load a 
helicopter could carry versus the gross weight of the 
helicopter. Fig. 1 shows data for thirteen different full 
scale helicopters plus the Bruiser. 

A linear trend analysis produced Eqn. 1: 



Payload = 27 0 + gross Wei ^ ht 

J 4.52 



( 1 ) 



The slope of the line was assumed to be true for all weight 
ranges, but the y- intercept was decreased for gross weight 
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Load vs Gross Weight 
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Figure 1 Payload vs Gross Weight 



less than 500 pounds. The final equation that was 



determined is shown in Eqn. 2. 



Gross Weight = 4 . 52 * Payload (2) 




The needs as stated in the previous chapter required a 
payload of between 20 and 30 pounds. This equates to a 
helicopter between 90 and 135 pounds. 
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The next trend that was studied was takeoff weight 
versus usable power. The desire was to use this trend to 
provide a general idea of what power would be required based 
on the above weight range. The trend shown in Fig. 2 was 
used to derive Eqn. 3. 

Usable Power = -70.66 ♦ 

7.38 



Takeoff Wt vs. Usable Power 
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Figure 2 Usable Power vs Takeoff Weight 
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Again, assuming the slope was accurate over the entire 
weight range, the intercept point was dropped to zero. The 
final equation is shown below. This equation shows that 
power required would range between 12 to 18 horsepower. 

usable Power - Take ° ff Wel ^ ht (4) 

7.38 

Next, the main rotor blade solidity and blade loading (Ct/cr) 
versus takeoff weight were compared. Fig. 3 shows that main 
rotor blade solidity is independent of takeoff weight. The 
minimum solidity was found to be 0.03 for the Robinson R-22 
Beta. The maximum solidity was 0.098 for the Aerospatiale 
AS 332 LI. The average solidity was approximately 0.07. 

Fig. 4 shows a similar trend to that of the blade solidity. 
The blade loading was also independent of takeoff weight. 
The minimum was 0.05 for the Bell 412. The maximum was 
0.095 for 

the Aerospatiale SA 365N, and the average was 0.075. 

The typical disk loading (DL) for model helicopters is 
between 1.0 and 2.0. Fig. 5 shows that the model helicopter 
does not fall within the normal range of full scale 
helicopters. The minimum full scale DL of 2.8 is that of 
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Takeoff Wt vs. Blade Solidity 
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Figure 3 Solidity vs Takeoff Weight 
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Takeoff Wt vs. Blade Loading 
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Figure 4 Blade Loading vs Takeoff Weight 
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Takeoff Wt vs. Disk Loading 
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Figure 5 Disk Loading vs Takeoff Weight 




the Robinson R-22. Full scale trends were inconclusive in 
providing guidelines for the model helicopter. Fig. 6 shows 
how power loading (PL) varies with disk loading. Power 
loading tells how much weight can be lifted for a given 
horsepower. The lower the DL, the greater the PL that can 
be achieved. In summary, the trend analysis shows that a 
desired payload of 20 to 30 pounds will require a helicopter 
weight between 90 and 135 pounds. An engine of 12 to 18 
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Disk Loading vs Power Loading 
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Figure 6 Power Loading vs Disk Loading 

horsepower would be required for a helicopter in this weight 
range. The solidity should be in the vicinity of 0.07 and 
the BL should be in the vicinity of 0.075. These trends 
provide a good starting point for the design process and a 
very good point of comparison for the preliminary design 
results . 
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IV. 



PRELIMINARY DESIGN 



A. CONSTRAINT DIAGRAM 

The first task to accomplish was to construct a rotor 
blade tip speed constraint diagram. Based on industry- 
criteria [Ref. 2:p. 90], a tip speed constraint diagram was 
made using the following limits. The upper boundary can 
either be set by a noise limit which is 750 fps or a hover 
tip mach number limit which is M<0.69 (771 fps). The hover 
tip mach number was used for this constraint diagram. The 
lower limit is set by requirement to store kinetic energy in 
the rotor system in case of power failure, in other words an 
autorotational limit. This limit is set at 400 fps. Two 
other limits that are important to consider are that of the 
advancing blade tip mach number limit and the advance ratio 
ifi) limit. The advancing blade tip mach number limit is set 
at M<0.8 to avoid compressibility effects on the advancing 
blade in forward flight. The advance ratio is set at fi< 0.4 
to avoid retreating blade stall at maximum forward speeds. 
One additional limit was incorporated to show the lower 
Reynolds number limit of 1.5 million. Fig. 7 depicts the 
constraint diagram developed using the above parameters. 
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Tip Speed Constraint Diagram 

Hummingbird 




Hover Limit M < 0.69 
Fwd Fit Limit M < 0.8 
Adv Ratio Limit Mu < 0.4 



Figure 7 Main Rotor Blade Tip Velocity Constraint Diagram 



B . PRELIMINARY ESTIMATES 

1. Gross Weight and Blade Radius 

Ref. 3 outlines basic steps to quickly estimate the 
gross weight and rotor radius. Using these steps, the 
following calculations were made: 
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Fuel Required 

(0.5 lb/hp-hr) Piston Engine 
Mission time 60 minutes 

Fuel required = (1.0) (25) (0.5)= 12.5 lbs of fuel* 
*This is approximately two gallons. 

Usable Load (UL) 

UL = crew + payload + fuel 
UL = 0 + 30 + 12.5 = 42.5 lbs 



Assuming 1970 technology, Fig. 8 shows a useful load 
per gross weight factor of 0.4 [Ref. 3:p. 641]. 

Gross Weight = 90.6 lbs. 
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Figure 8 Historic Trends of Ratio of Useful Load to Gross 
Weight 



This number for gross weight falls within the range 
developed by the trends analysis, though it is on the low 
side for a payload of 30 pounds. Realizing that the design 
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weight of any design will grow, it is assumed that the 
higher weight is more reasonable; therefore, a 130 lb gross 
weight will be used for the following estimates. Based on 
this weight and the DL range from 1.0 to 2.0, the rotor 
radius ranges from 4.55 to 6.43 feet. This rotor radius 
range is validated by calculation of required blade area as 
described in Ref. 2. Assuming a maximum forward velocity 
( V m ax ) knots, Fig. 9 shows that blade area required is 

5.5 square feet, which is 1.83 ft 2 per blade for a three 
blade system and 1.375 ft 2 for a four blade system. This 
translates into a radius of 4.9 ft, a DL of 1.74 lb/ft 2 with 
a chord of 0.375 ft (4.45 in) for the three blade system, 
and a radius of 4.7 ft, a DL of 1.86 and a chord of 0.292 ft 
(3.5 inches) for the four blade system. 



■t* V 




Figure 9 Determination of Blade Area for New Rotor Design 
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Either of the above systems would be feasible for this 
design, but the four blade system does have two distinct 
advantages. First, the blade vibrations in a four blade 
system are less than those in a three blade system. Second, 
the aspect ratio (A R) of the three blade system is 13.1 and 
the AR for the four blade system is 16.1. The normal range 
for AR is 15 to 20. Any blade below that will have a lower 
blade efficiency due to tip losses and any blade above that 
could pose structural problems. 

2 . One Hour Estimate 

In this estimate, there are two constraints: hover 
out of ground effect (HOGE) and high forward speed flight. 
The following is a list of the tentative performance 
requirements : 

Payload 20-30 lbs 

Crew 0 lbs 

Max Speed @ S.L. 70 knots 

Cruise Speed (0 . 9*V max ) 63 knots 

Vertical Rate of Climb 450 fpm @ 4000 ft 95 deg. F 
Engine: One - Max continuous rating 25 BHP. 

The following calculation as detailed in Ref. 4 will be done 
using the previous gross weight estimation of 130 lbs. Fig. 
10 shows that for a weight of 130 lbs, the equivalent flat 
plate area (f e ) would be approximately 0.8 ft 2 . [Ref. 5:p. 

35] 



18 




Figure 10 Parasite Drag Area vs Design Gross Weight 

a. Maximum Forward Velocity ( V~ mi 1 f ) 

Assuming seventy percent of installed power is 
used to overcome parasite drag at high speeds, maximum sea 
level speed can be calculated. 

_ r (30 min rating eng) •, j (5) 

’'max * x 1 f 1 

A e 

Maximum BHP = 25 

Assuming Figure of Merit of 70% 

THP = 17.5 
V v = 114 knots 

max 

Based on the assumption that the helicopter was to be 
derated to a power of 18 HP, the calculations would be 
changed as follows: 
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Maximum BHP 
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Assuming Figure of Merit of 70% 

THP = 12.6 

V max = 41* [12. 6/. 8] A (1/3) 

V = 102 knots 

max 

b. Rotor Sizing 

The preliminary size of the rotor is determined 
using the rate of climb (ROC) stated earlier. Additional 
power required to climb 450 fpm is approximately 10% of the 
power required for HOGE. 

ROC 450 fpm @ 4000 ft 95 deg. F (95% of 17.5 THP) 

95% of 17.5 hp = 16.63 hp 

PL = 130 lb/16.63 hp = 7.82 

From Fig. 11 the DL = 3.5 [Ref. 4:p. 8]. 

DL = Weight/ (pi*R A 2) =3.5 
R = 3.44 ft 

Again assuming a derated engine, 

95% of 12.6 hp = 12.0 hp 
PL = 130 lb/12.0 hp = 10.83 
From Fig. 11 
DL = 1.75. 

DL = Weight/ (pi*R x 2) = 1.75 
R = 4.86 ft 
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Figure 11 Rotor Performance for Design Conditions 
c. Starting Point 

The preceding calculations provide a gogd point 
at which to begin more in-depth design. In summation, the 
starting point is 

Gross Weight 130 lbs 
Empty Weight 87.5 lbs 
Payload 30.0 lbs 

Fuel Weight 12.5 lbs 
Rotor Radius (R) = 4.86 
Disk Loading (DL) = 1.75 
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3 . Design from Scratch or Procure 

Based on this data, there were two choices to make. 
First, design a helicopter from scratch which would meet the 
needs of the Aeronautics and Astronautics Department. While 
this is definitely an exceptional learning process, the end 
result would leave the department with a paper helicopter 
and at least a year's worth of manufacturing. The other 
path of choice was to look for a helicopter on the open 
market that closely approached the stated needs. The second 
alternative, being the most productive and time smart, was 
chosen. Though the possibility of locating a helicopter 
that would meet every need was slim, it was the best choice, 
realizing that the redesign of an existing, flying 
helicopter was much more frugal financially and time-wise. 
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V. 



HUMMINGBIRD ACQUISITION 



A. CAPABILITIES 

A helicopter of the correct size and payload capability 
was located at the GMP model helicopter company. The owner 
and designer, Mr. John Gorham, the lead engineer with 
Lockheed Corporation on the L-1011, had built ten 165 pound 
RPVs for the U.S. Army. Their original purpose was to be 
used as Soviet Hind-D recognition devices. This vehicle was 
nearly ideal for meeting the needs of the department. There 
were some shortfalls, which will be mentioned in the next 
section, but nothing of a critical nature that could not be 
redesigned or changed. The department purchased one of Mr. 
Gorham's helicopters, with the desire to purchase a second. 
The initial RPH was named Hummingbird I. The purpose of 
buying two helicopters is to allow for concurrent HHC, 

NOTAR, and ALTOS research. One helicopter would be used in 
its original configuration for NOTAR and ALTOS research. 

The other would be converted into a platform capable of HHC 
research. Fig. 12 contains two pictures of the Hummingbird. 
Table I contains a list of the Hummingbird characteristics 
and capabilities as it was received from Mr. Gorham. 
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Table I HUMMINGBIRD CHARACTERISTICS AND CAPABILITIES 



Hummingbird : 




Characteristics 




Weights 




Max Gross Weight 


165 lbs 


Empty Weight 


115 lbs 


Fuel capacity 


6.5 lbs 


Rotor Parameters 


Main Tail 


Radius (R) 


5' 3" 12.5" 


Chord (c) 


6" 2.625" 


Solidity (sigma) 


! 0.0061 0.02 


No. of blades (b) 


2 3 


Tip speed 


! 303 fps 241 fps 


Twist 


: -5 0 


Hinge offset ration (e/R) 


0.127 0.24 


Airfoil 


NACA 0012 NACA 0012 


Engines 




Type 


Westlake 342 Series 2100D 


Number 


1 


Maximum Usable Power 


25 BHP @ 7000 rpm 


Maximum Torque 


25 ftlb @ 4000 rpm 



B . SHORTFALLS 

The Hummingbird satisfied many of the department's 
needs; specifically, it was large enough to be developed 
into a quarter scale model, it had a 50 plus pound payload, 
and in addition, it had already proven itself in flight 
test. There were, however, some shortcomings that had to be 
addressed. These shortfalls included the need for more than 
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a two blade hub, an RPV that possessed autorotational 
capabilities and could produce blade Reynolds numbers which 
would be comparable to full scale helicopters. 

1. Two Blade Rotor 

In order to have the capability to do HHC research 
work on the RPV, it was required that the helicopter have at 
least three rotor blades or more. The reason for this 
requisite is that the vibratory forces at the rotor blade 
root are produced by the (n-1) /rev, n/rev, and (n+1) /rev 
vibrations. The n/rev vertical forces and moments are 
transmitted to the fixed system at a frequency of n/rev. 

The (n-1) /rev and (n+1) /rev flapwise blade root shears 
result in n/rev hub pitching and rolling moments in the 
airframe. The n/rev flapwise blade root shears feed into 
the airframe as n/rev vertical forces. The (n-1) /rev and 
(n+1) /rev chordwise root shears produce n/rev airframe hub 
forces in the fore and aft and lateral directions. The 
n/rev chordwise root shears result in n/rev hub yawing 
moments. Choosing a rotor system with three or more blades 
therefore will not interfere with the primary 1/rev control 
inputs. [Ref. 6:p 19] This change alone created a 
significant amount of work because it includes changing the 
main rotor hub and redesigning the rotor blades. 
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2. Lack of Autorotational Capability 

The original design never incorporated the ability 
for the Hummingbird to autorotate, and this was considered 
unacceptable for the NPS flight research vehicle. The 
bearings in the drive system do not allow the main rotor to 
freewheel should the RPV have an engine failure. Presently, 
should the engine fail it will cause a rapid decay in rotor 
rpm which will transform the RPV from a flying machine to a 
projectile. This is a critical redesign requirement which 
must be accomplished in order to protect the Navy's 
financial and research investment in the Hummingbird. 

3 . Low Reynolds Number 

In the HHC research, it is not enough to prove that 
the concept will reduce vibration; that result is 
sufficiently proved. The desire is to have the ability to 
test whether HHC also provides performance improvements. 

This comparison can only be accomplished on helicopters with 
similar Reynolds numbers. 

Reynolds number is the ratio of the inertia to the 
viscous forces on a volume of fluid (Eqn. 6) . 

Re = A™. (6) 
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The importance of Reynolds number in the comparison of 
geometrically similar bodies in incompressible flows can be 
shown using the Navier-Stokes equations. Eqn. 7 shows the 
dimensionless Navier-Stokes equation. 

i!iL'=-.§£-'._Lv'V- 

Dt‘ Sx' Re Dt' ay' Re Dt< 3z' Re O) 

These equation show that, given geometrically similar 
bodies, the equations of motion are identical for the same 
Reynolds number. The similar bodies includes surface 
roughness as well as shape. [Ref. 7:p. 304] 

The Hummingbird's Reynolds number was on the order of 
0.9 million, and it needed to be in the range of 1.5 to 1.8 
million to be comparable to that of the 0H-6A's Reynolds 
number equal to 2.4 million. 
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VI. 



BLADE DESIGN 



A. SPECIFICS 

1. Airfoil 

The airfoils that were available with the 
Hummingbird were NACA 0012 and NACA 0013 . This type of 
airfoil is widely used in industry for numerous reasons. 
First, the airfoil is symmetrical and therefore there is no 
nose -down pitching moment which is associated with cambered 
airfoils. Second, it is a relatively thick airfoil which 
will provide an acceptable maximum lift coefficient. 

Finally, there is a vast amount of data available on this 
airfoil, which provides for much easier analysis. This is 
due to the fact that the NACA 0012 airfoil was selected by 
Mr. Sikorsky for the VS- 300, the world's first truly 
successful helicopter. It was also the airfoil of choice 
for almost all early helicopters, including more recent 
aircraft still in service, such as the Navy H-3. 

For improved performance, a more recent advanced 
technology airfoil worthy of consideration is the NACA 
23012. The characteristic drooped nose is an effective 
method of increasing the maximum lift coefficient. Also, at 
high lift it has lower drag than a similar six- series 
airfoil at low mach numbers. [Ref. 3:p. 388] 
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2. Twist 



The twist of a rotor blade enhances two main areas 
of the helicopter's performance- -hover performance and 
retreating blade stall- -but it also produces increased 
vibration in forward flight. Its most notable adverse 
effect is to increase blade vibratory stresses and in this 
way reduce blade fatigue life. Built-in blade twist affects 
the radial variation of inflow angle from blade tip to blade 
root. Ideal twist is represented by Egn. 8. 

* - *« B <8) 

r 

Hovering performance is enhanced with the addition of 
negative twist by creating a more uniform inflow 
distribution along the blade span. The larger the amount of 
twist, the closer it approximates an ideal twist 
distribution. Generally accepted values of main rotor blade 
twist are -8 to -14 degrees. Ideal twist yields the minimum 
induced loss for a given thrust. Retreating blade stall is 
delayed when twist is employed by unloading the tips, which 
reduces the tip angle of attack [Ref. 8:p. 57]. 

These regions of enhanced performance are not 
critical areas for the Hummingbird. The Hummingbird will 
operate at relatively slow forward speeds; therefore, 
retreating blade stall is not a concern. It would be 
attractive to have highly efficient blades in the hover 
regime, but with the present payload capability of 50 
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pounds, it is not critical. If these were critical areas 
for the Hummingbird, it might be worth the effort to 
manufacture new blades with twist; but since they are not, 
simplicity rules, and it was decided not to incorporate 
twist in the present rotor blade design. 

3. Tip Velocity 

Determining the tip speed (V tj ) for the Hummingbird 
was a difficult matter because its effects were coupled with 
so many other areas of the helicopter's performance. The 
constraint diagram from Chapter II is shown again in Fig 13 . 

The design point is shown to be in the proximity of V max of 
40 knots and V. . of 450 f ps . V mav was obtained from the 
Hummingbird original design. The means by which V tip was 
determined will follow. Low tip speeds have the advantage 
of low noise and good hovering performance. "High tip 
speeds have the advantage of low rotor and drive system 
weights and high stored energy for autorotative entries and 
flares." [Ref. 3] In the case of the Hummingbird, V tjp was 
the variable used for a tradeoff study. V tj had to be large 
enough to produce Reynolds numbers which could be compared 
with full scale helicopters, but small enough to provide a 
reasonable Figure of Merit (FM) . 
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Tip Speed Constraint Diagram 

Hummingbird 



1,000 




Hover Limit M < 0.65 
Fwd Fit Limit M < 0.6 
Adv Ratio Limit Mu < 0 4 



Figure 13 Constrain Diagram 



B . TRADEOFF 

The need to achieve a Reynolds number in the range 
mentioned earlier was the source of some consternation. If 
Reynolds numbers were not considered, it was quite easy to 
design a rotor blade that would fly efficiently and meet the 
department needs. The complication begins when the V tjp , 
Reynolds number, Figure of Merit and chord length 
requirements are all met at the same time. Reynolds number 



32 



is a function of V u and chord length as shown by Eqn. 9. 

Re = 6400 V tip c ( 9 ) 

Fig. 14 and 15 show this relationship over a certain range 




of V Up and chord length. 

Figure of Merit is shown by the relationship in Eqn. 10. 
Induced power is the power required to overcome the drag due 
to the generated lift. Total power also includes the power 
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FM = 



( 10 ) 



P_i 

D 

total 

required to move the rotor blades through the air (profile 
power) and the power required to drag the fuselage through 
the air (parasite power) . Fig. 16 and 17 show this 
relationship over the same V tj and chord length range. 

To link the Figure of Merit and Reynolds number together, 
lines of constant Reynolds numbers were plotted on the 
Figure of Merit versus chord length chart. Fig. 18 allowed 
for the tradeoff to be visualized simultaneously on one 
graph. The Hummingbird came with three sets of blades. One 
set was of radius 4.52 feet with a chord of 6 inches. The 
other two sets were of radius 5.0 feet with a chord of 6.5 
inches. To expedite the design process, the available 
blades were selected versus designing and manufacturing new 
ones. Since there was a requirement for at least a three 
bladed rotor, the second set of blades was chosen. 

The advantages of this choice are seen in the tradeoff 
study. Choosing a longer chord length allowed for a slower 
tip speed for the same Reynolds number. Since profile power 
is proportional to V tfp 3 , the slower tip speed increases the 
FM from 0.48 to 0.54. 
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C. STABILITY 

The pitch and roll damping is produced by the tilt of 
the tip path plane, which lags behind the motion of the 
shaft by an amount that is proportional to the rate of pitch 
or roll. Therefore, the aerodynamics on the blade causes 
the tip path plane to tend to stabilize itself in an 
equilibrium position with respect to the shaft. "If the 
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aerodynamic and inertia flapping moment are equated, the 
following results for the angular displacement of the rotor 
plane with respect to the shaft per unit tilting velocity of 
the shaft is obtained for the hovering case: "[Ref. 8:p. 275] 




A = JJL 

oo yCl 



( 11 ) 



The quantity 16/ (yfl) can be understood more easily by 
examining its physical interpretation as follows: The thrust 
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vector lags the rotor shaft by a time constant of 16/(70) 
seconds if the rotor shaft is tilted with any constant 
angular velocity. Therefore, the larger the time constant 
(16/(70)), the greater the system is damped. Lock number 



Tradeoff Study 

Chord vs Figure of Merit 




4.0 4.5 5.0 5.5 6.0 6.5 

Chord (in) 

\ 

Blade Tip Velocity (fps) 

— 300 400 *450 -*500 * 550 * 600 

V ) 

Bladed Rotor System 



Figure 17 Figure of Merit vs Chord Length 

(7) relates the inertia and aerodynamic characteristics of a 
blade and is shown in Eqn. 12. 

Y = C P (12) 
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Figure 18 Figure of Merit and Reynolds Number Tradeoff 



The stability of the Hummingbird is a critical area that 
will require additional research. The two proposed 
modifications to the helicopter significantly alter the 
stability of the RPV. Eqn. 11 shows that rotor speed and 
Lock number will affect the damping of the rotor. The rotor 
blade is not changed; therefore, the Lock number will remain 
the same. The rotor speed will be increased from 550 rpm to 
759 rpm, thus decreasing the effective rotor damping. Ref. 

8 states that 
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In addition to the effects of rotor speed, rotor damping 
may be increased by the use of devices that act upon the 
control system in such a manner as to increase the 
displacement of the rotor from its trim position due to 
a given rate of roll or pitch. An example of such a 
device is a rate gyro that would apply opposite control 
by an amount proportional to the rolling or pitching 
velocity of the helicopter. 

The Bell -Hiller stabilizer bar hub configuration is such a 
device. Therefore, changing from that type of hub to an 
articulated hub will significantly decrease the rotor roll 
and pitch damping and thus the helicopter's handling 
characteristics. For the 3-bladed Hummingbird, rate gyros 
may be required in order to have satisfactory flying 
qualities . 

D. BLADE DYNAMICS 

1. Introduction 

It is very important to consider the vibratory 
effect on the main rotor blades. High vibrations result in 
high vibratory shear at the rotor hub, which results in high 
fuselage vibrations. The high vibrations also cause 
vibratory shear stress, which reduces the effective life of 
a rotor blade. The vibratory resonance response also 
affects dynamic stability of the rotor system. 



39 



2 . Main Rotor Blades 

The following analysis was done on two different 
rotor blades. The first was termed the "heavy" blade. It 
had a length of 4 ft 6.25 inches, chord of 6 inches, 
thickness ratio of 12 percent and a weight of 5.18 lbs. The 
second blade, termed the "light" blade, has a length of 5 
ft, chord of 6.5 inches, thickness ratio of 13 percent and a 
weight of 4.66 lbs. The length given for the blades above 
does not include the 8 inch blade offset of the rotor hub, 
but it was included in the rotor radius when the analysis 
was performed. 

In order to analyze each blade, they were broken up 
into segments: the heavy blade into 17 segments and the 
light blade into 24. Fig. 19 and 20 show a schematic for 
the two blades and how the blades were divided into 
segments. The analysis for each blade was conducted with 
the exact same technique, and therefore the discussion will 
be limited to the heavy blade. The following quantities 
were required to perform the blade analysis: section radii, 
chord, segment width, segment volume, segment weight, and 
area moment of inertia (Ixx) , and material modulus of 
elasticity. The section radii was the distance from the 
center of rotation to the center point of the segment. The 
chord was constant for the blade until the small taper at 
the root. The segment width was constant with the exception 



40 



of the tip segment, which was allotted any leftover blade. 
Fig. 21 shows a cross section of the blade, which was used 
for area and moment of inertia calculation. For 
calculations, the airfoil cross section was broken up into 
three sections and approximated as follows: (1) the nose 

section by a parabola; (2) the center section by a box; and 
(3) the tail section by an isosceles triangle. 

The cross sectional area of each section was 
calculated using standard geometric formulas in Eqn. 13. 

Parabola: A = Box : A = bh Triangle: A = bh 

3 2 

The area moments of inertia for the box and triangle were 




41 




calculated using known formulas, while integration was used 
to determine the Ixx of the parabola. (Eqn. 14.) 

Box : Ixx =-^-bh i Triangle-. Ixx = -^-bh* Parabola : Ixx = f y 2 dA 
12 3 12 Ja 

( 14 ) 

The parabola was approximated by the equation x = 6.67y 2 . 

Using the calculated cross sectional area and the 
segment width, the volume was calculated. The segment 
volume was then used to determine the segment weight. The 
assumption was made that the rotor blade was made of a 
homogeneous material - spruce. The only variation in these 
calculations was where the blade was tapered very near the 
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root. This section was no longer an airfoil section, but an 
inch thick tapered block. The segment weight and Ixx of the 
root section also included contributions from the blade grip 
assembly. Table II shows the data for each blade. 
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Table II BLADE STATION DATA 



Heavy Blade Radius « 4 52 It 

Weight 5 18 lbs Chord = 6 in 
t/c = 12% 



Section 


Radii 


Chord 


Xaec Ar Sea W4 Volume Seg Wt 


Ixx 


Pelt* Cf 




(ini 


M 


{ialfl 


liQl 


liar. 3) 


(M 


Halil 


mibi 


1 


59 1 


80 


1 68 


6 25 


1050 


0459 


0 097 


230 85 


2 


54 5 


60 


1 68 


300 


504 


0 220 


0 097 


49 03 


3 


51 5 


60 


1 68 


300 


5 04 


0 220 


0 097 


46 33 


4 


48 5 


60 


1 68 


300 


504 


0220 


0 097 


43 63 


5 


45 5 


60 


1 68 


300 


504 


0 220 


0 097 


40 93 


6 


425 


60 


1 68 


300 


5 04 


0 220 


0 097 


38 23 


7 


39 5 


60 


1 68 


300 


504 


0 220 


0 097 


35 53 


8 


36 5 


60 


1 68 


300 


504 


0 220 


0 097 


3 Z 64 


9 


33 5 


60 


1 68 


300 


504 


0 220 


0 097 


30 14 


10 


305 


60 


1 68 


300 


5 04 


0 220 


0 097 


27.44 


11 


27 5 


60 


1 68 


300 


504 


0 220 


0 097 


24 74 


12 


24 5 


60 


1 68 


300 


5 04 


0 220 


0 097 


22 04 


13 


21 5 


60 


1 68 


300 


504 


0 220 


0 097 


1934 


14 


18 5 


60 


1 68 


300 


5 04 


0 220 


0 097 


1664 


15 


155 


54 


5 40 


300 


1620 


0 708 


0 450 


44 82 


16 


125 


50 


500 


300 


1500 


0 656 


0417 


33 47 


17 


95 


38 


3 75 


300 


11 25 


0 742 


0 608 


28 77 



Light Blade Radius « 5 0 ft 

Weight = 4 66 lbs Chord = 6 5 in 
t/c = 13% 



Section 


Radii 


Chord 


Xaec Ar Seg Wd Volume Seg W|. 


Ixx 


Delta Ct 




m 


(ml 


(m: a 


(in) 


tiaia 


m 


Mill 


(ftlbl 


i 


65 6 


65 


3 37 


48 


16 18 


0 353 


0 478 


151 24 


2 


62 0 


65 


3 37 


24 


609 


0 177 


0476 


3574 


3 


59 6 


65 


3 37 


24 


809 


0 177 


0478 


34 35 


4 


57 2 


65 


3 37 


24 


809 


0 177 


0478 


3297 


5 


54 8 


65 


3 37 


24 


809 


0 177 


0 478 


31 59 


6 


52 4 


65 


3 37 


24 


8 09 


0 177 


0 478 


30 20 


7 


500 


65 


3 37 


24 


809 


0 177 


0478 


28 82 


8 


47 6 


65 


3 37 


24 


809 


0 177 


0478 


27.44 


9 


45 2 


65 


337 


24 


809 


0 177 


0 478 


26 05 


10 


42 8 


65 


3 37 


24 


809 


0 177 


0478 


24 67 


11 


40 4 


65 


3 37 


24 


809 


0 177 


0 478 


23 29 


12 


380 


65 


3 37 


24 


809 


0 177 


0 478 


21 90 


13 


35 § 


65 


3 37 


2 4 


809 


0 177 


0478 


2052 


14 


332 


65 


3 37 


24 


809 


0 177 


0 478 


19 14 


15 


308 


65 


3 37 


24 


809 


0 177 


0.478 


17.75 


16 


28 4 


65 


3 37 


2 4 


809 


0 177 


0478 


1637 


17 


26 0 


65 


337 


24 


809 


0 177 


0478 


14 99 


18 


23 6 


65 


3 37 


24 


609 


0 177 


0478 


1360 


19 


21 2 


65 


3 37 


2 4 


809 


0 177 


0478 


1222 


20 


188 


65 


3 37 


24 


809 


0 177 


0 478 


1084 


21 


164 


65 


3 37 


24 


809 


0 177 


0478 


9 45 


22 


140 


58 


5 76 


24 


13 82 


0 302 


0 480 


13 79 


23 


11 6 


48 


480 


24 


11 52 


0 252 


0 400 


952 


24 


92 


42 


420 


24 


1008 


0 470 


0 845 


1411 
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3 . Myklestad Determinant Method 

There are many methods for analyzing vibrations in a 
rotor blade. Some of these are 

• Rayleigh-Ritz 

• Holzer 

• Myklestad- Prohl 

• Matrix iteration 

The Myklestad method, which determines the natural 
frequencies and modes of the blades, was the method of 
choice. It is especially well suited for vibration analysis 
of rotor and turbine blades. Dr. Nils Myklestad was a 
consultant to Bell Helicopter Company for many years. The 
proper analysis requires considering the coupled flapwise- 
edgewise- torsional response of the blades. Because the 
edgewise and torsional stiffness is much greater than the 
flapwise stiffness and thus has much high frequencies, only 
the flapwise responses of the blades were considered. 

The uncoupled flapwise Myklestad system is shown in 

Fig. 22. 

From this diagram the equations of equilibrium for this n th 
element can be written. 
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Centrifugal Force : 



Tn - T m+1 +11^ 2 r N 



( 15 ) 



Shear: 



S N =S N + i +m N <*> 2 Yn 



( 16 ) 



Moment : 



Mn = 1 n.n + 1-'WyN + 1-y N ) 



( 17 ) 
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Slope : 



9 N = e N+1 d + T N+1 i^)-M N+1 ^2-s 



N+1 * 



L n,n+1 

2EI 



(18) 



Deflection: 



Yn = y»M- 0 N 1 n,n*i + 



^N+1 ®N+1 



n,n+1_^N+1^- n,n+1_^N+1^- n,n+1 



(19) 



3EI 



2EI 



3EI 



The natural frequencies of the system are obtained 
by assuming a frequency (Q) and then calculating the 
centrifugal forces, shears, moments, slopes, and deflections 
of each blade segment from the blade tip to the blade root. 
The boundary conditions at the blade tip are 
S T = Mj = 0; 0 T = 0 T ; y T = y T 
At the root of the blade, the two unknowns 0 T and y T are 
carried along and the equilibrium equations are written as 
follows : 

S o = AsYt + B s e r 

M o = + b m 0 t 

9 o = + B 9 0 T 

y 0 = Vt + B y 0 T 

In the preceding equation, the coefficients A and B will be 
calculated as the process proceeds down the blade. They are 
functions of mass and stiffness properties of the blade, the 
rotational speed, and the assumed frequency. The boundary 
conditions at the root of the hinged blade are 
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0 



M o - y 0 = 
or 

+ B m 0 t = 0 
Vt + B y 0 r = 0 

These equations can be written in matrix form. The 
frequencies that satisfy the equation (when the determinant 
equals zero) are the natural frequencies and can readily be 
determined. 

A Fortran code, Appendix A, was written to solve the 
Myklestad equations and determine the flapwise natural 
frequencies of both Hummingbird rotor blades. Fig. 

23,24,25, and 26 show a graphical output of the Myklestad 
determinant for both blades. 

The natural frequencies, where the line crosses the 
x-axis, are clearly seen. Appendix B contains a numerical 
output of the same data for the different natural 
frequencies . 

A means to verify the accuracy of the code's ability 
to determine the natural frequencies was to check the ratio 
of the determinant coefficients. Ay and A^ are the 
displacement and slope coefficients respectively for the 
rotor blade root section in response to starting boundary 
conditions at the blade tip of slope equal zero and 
deflection equal one. B y and B m are similar to Ay and A^ 
except the starting boundary conditions are slope equals one 
and displacement equals zero. Knowing that the determinant 
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Myklestad Determinant 

Hummingbird 




Operating rpm 759 
Heavy Blade 



Figure 23 Myklestad Determinant (Heavy Blade) 
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Myklestad Determinant 

Hummingbird 



20 



15 



o) 

c 



~ 5 

d) 

Q 



-5 




100 200 300 400 500 

Omega (rad/s) 



600 



700 



Operating rpm 759 
Heavy Blade 



800 



Figure 24 Myklestad Determinant (Heavy Blade) 
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Myklestad Determinant 

Hummingbird 




Operating rpm 759 
Light Blade 



Figure 25 Myklestad Determinant (Light Blade) 
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Myklestad Determinant 

Hummingbird 




Operating rpm 759 
Light Blade 

Figure 26 Myklestad Determinant (Light Blade) 

equals zero at the natural frequency, the ratio of -A y /B y 
should equal -A^B^ Using data in Appendix B for the heavy 
blade's first natural frequency, the accuracy is seen. 

1st Natural freq. = 107.1522 rad/s 
Ay= 7.48560 B y = -469.659 

A^= 1 . 79121E+04 B m = -1.12384E+06 

-Ay/By = 1 . 593837E- 02 - A m /B |n = 1.593830E-02 



52 



It is interesting to note that the natural 
frequencies are much higher than those of conventional 
helicopters. This is attributable to the fact that model 
rotor blades are much shorter and made of solid spruce, 
which makes them considerably stiffer than full scale 
blades . 

4 . Blade Modes 

The blade modes were used to verify that the 
calculated frequencies were in actuality the ascending 
natural frequencies. Appendix C contains a Fortran code 
written by Lt . M. Avila that was modified for use on the 
Hummingbird rotor blades to obtain the rigid and bending 
flapwise mode shapes. Fig. 27 shows the mode shapes for the 
heavy blade while Fig. 28 shows the mode shapes for the 
light blade. 
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Rotor Blade Analysis 

Hummingbird 




9.5 14.5 19.5 24.5 29.5 34.5 39.5 44.5 49.5 54.5 

Blade Station (inches) 

/• \ 

Bending Modes 

Rigid Flapwise — 1st Flapwise ^2nd Flapwise "5- 3rd Flapwise 

V ^ 

Heavy Blade 

Figure 27 Heavy Blade Mode Shape 



5. Southwell Plot 

The Southwell plot, also called a Fan plot, is a 
resonance diagram of the individual blade. It is used to 
depict the flapwise modes in relation to the integer 
multiples of the rotor rpm, commonly referred to as IP, 2P, 
3P etc. For a good design, effort should be directed toward 
keeping the blade natural frequencies away from the 1P,2P, 
etc. forcing frequencies, since these are the primary 
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Rotor Blade Analysis 

Hummingbird 




9.2 14.2 19.2 24.2 29.2 34.2 39.2 44.2 49.2 54.2 59.2 64.2 



Blade Station (inches) 



Bending Modes 

Rigid Flapwise “1st Flapwise ^ 2nd Flapwise -5* 3rd Flapwise 

\ J 

Light Blade 

Figure 28 Light Blade Mode Shape 

excitations. The Fortran code used for the Myklestad 
determinant was run for rotor speeds ranging from 0 to 90 
rad/s. The results of these runs are shown in Table III. 
The area of concern is within plus or minus 6 rpm of the 
operating rpm. Fig. 29 and 30 show the Southwell plot for 
the heavy blade and light blade respectively. 

The rigid mode always lies very near the IP line. 
The flapwise bending modes for both blades are much higher 
than that for full scale blades. Again, this is attributed 
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Table III BLADE FLAPWISE RESPONSE 



Heavy Blade Light 

Blade 


rad/s 


rpm 


Rigid 


1st 

Bend 


2nd 

Bend 


Rigid 


1st 

Bend 






(cpm) 


(cpm) 


(cpm) 


(cpm) 


(cpm) 


0 


0 


1.29 


1758.9 


5757.3 


1 . 19 


2876.2 


10 


95.5 


129.9 


1791.5 


5784.0 


119 . 7 


2893.8 


20 


190 


262.6 


1887.3 


5863.3 


240.6 


2945.7 


30 


286.5 


390.6 


2036.3 


5992.7 


360.7 


3030.0 


30 


382 


517.6 


227.9 


6169.8 


480.9 


3144.2 


50 


477.5 


643.6 


2452.3 


6389.4 


601.1 


3285.4 


60 


57.3 


774.4 


2700.1 


6647.3 


721.9 


3449.2 


70 


668.5 


902.4 


2967.4 


6939.5 


840.9 


3634.5 


79.42 


759 


1023 . 7 


3229.2 


7241.9 


953 . 7 


3824.5 


90 


859.4 


1159.3 


3535.2 


7608.9 


1080 


4053.7 



to the greatly increased stiffness. The second and third 
inodes for the light blade were higher than those for the 
heavy blade. This is accounted for by the fact that 
although the light blade is larger, the lower mass more than 
offsets the length by increasing the natural frequency. The 
heavy blade Southwell plot show that there will be no 
vibratory resonance problems with any of the modes. The 
first flapwise bending mode of the light blade lies very 
near the 5P line at the operating frequency. The 5P lies at 
62.56 Hz while the second natural frequency lies at 63.74 
Hz. The close proximity of the two frequencies should not 
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be a problem for the Hummingbird, as the critical 
frequencies occur at the 2P, 3P, and 4P for a three blade 
rotor system. For a full-scale equivalent 3-bladed 
helicopter, this near 5P resonance might need to be 
monitored due to possible high amplification of 5P flapwise 
blade stresses. In the case of the Hummingbird, the 
relatively short blades are sufficiently overdesigned to 
avoid this problem. 

Discussion of blade resonance results for the RPH blades 
would not be complete without comparing these values with 
values obtained for full-scale helicopters. Here for the 
light and heavy blades, we find the first rigid mode at 
1.260 and 1.350 respectively. For a uniform articulated 
blade this would be considerably closer to IP at 1.020.- The 
first flapwise bending mode for an articulated fullscale 
blade would be in the vicinity of 2.50, whereas analysis 
shows in this case that we have much higher frequencies of 
5.03fi and 4.250 for the light and heavy blades respectively. 
Similarly, the second flapwise bending mode for a full-scale 
articulated blade is about 4.70 to 4.90, whereas for the 
heavy blade, the second flapwise mode occurred much higher, 
at 9.540. This is even higher than where we would normally 
expect the frequency of the third flapwise mode to occur for 
a full-scale rotor. Normally the third mode resonance 
occurs in the vicinity of 8P. 
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Figure 29 Southwell Plot (Heavy Blade) 
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Figure 30 Southwell Plot (Light Blade) 
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VII. 



PERFORMANCE ANALYSIS 



A. POWER CALCULATIONS 
1 . Power Required 

Using the parameters established in the blade design 
section and the procedures outlined in Ref. 9, power 
estimations were calculated. There exist two regimes of 
operation that establish the power required limits. The 
first is hover out of ground effect; the second is high 
speed flight. Since the Hummingbird would not operate in 
the high speed flight regime, the HOGE was established at 
the limiting flight condition. Two hover conditions were 
evaluated: sea level standard day and 4000 ft MSL 95 degrees 
F. As alluded to earlier, there are three basic power 
calculations that must be accomplished. They are induced 
power (Pj) , the power required to produce lift; profile 
power (P 0 ) , the power require to push the blades through the 
air; and parasite power (P p ) , the power required to move the 
rest of the fuselage through the air. The following 
formulas were used to calculate these quantities. 



rn 1 . 5 

HP; 

v/2p A B 550 




( 20 ) 
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C do p oA V ci p 



P n = 

° 4400 



HP 



( 21 ) 



= 






550 



HP 



( 22 ) 



T = Thrust B = Tip Loss Factor b = Number of blades 
A = Rotor Disk Area V f = Forward Velocity 

f e = Equivalent flat plate area 
Using the above formulas, the total power required to hover 
out of ground effect was calculated. Induced power was 4.95 
HP, while profile power was 4.4 HP. The parasite power is 
equal to zero in the hover case. To calculate the power for 
the 4000 ft 95 degree F case requires a density correction. 
The power required to hover in ground effect (HIGE) was also 
calculated. The induced power ratio equation, based on the 
a ratio of the rotor height above the ground over the rotor 
diameter, is used for this calculation, and is shown below: 



P • ( 7~ GE) h 4 h 3 h 2 h 

=-0 . 127 6 ( — ) +0 . 7 0 80 ( — ) -1.4569 ( — ) +1 . 3432 ( — ) +0 . 51 

P i ( OGE) D D D D 

(23) 

The rotor height above the ground (h) equals 2.5 feet. 

Using the light blades for these calculations, the ratio of 
h/D equals 0.22. It can be shown that the induced power 
required to hover IGE is 3.7 HP. These calculations are 
solely for the main rotor system. The power required for 
the tailrotor will be included in later calculations. The 
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final flight regime that must be considered is forward 
flight. The profile power in forward flight is shown by 
Eqn. 24. 



This equation does not include high speed effects such as 
compressibility or retreating blade stall, but as stated 
previously, these will not be factors for the Hummingbird. 
Power required calculations for velocities ranging from 
hover to 100 knots are shown in Fig. 31. 




(1+4 . 3 |X 2 ) 



(24) 
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Figure 31 Power Required Profile 
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Tailrotor power is included in these calculations, assuming 

I 

that it is nine percent of main rotor in the hover 
condition, decreasing to three percent at mid-range 
velocities. Appendix D contains the spreadsheet calculation 
for all power required calculations. 

2. Power Available 

The power available curve, Fig. 32, was constructed 
from data obtained from Ref. 10 and 11. An assumption of 
near linearity over the entire RPM range was made. 



Power Available 

Hummingbird 

WAE Limited 342 




Figure 32 Power Available vs RPM 
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Fig. 32 will be used again when determining the derated 
power that will be use in conjunction with resizing of the 
main rotor drive gear for obtaining desired rotor RPM. 

B. RETREATING BLADE STALL 

A phenomena that is particular to helicopters is 
retreating blade stall. It occurs at high forward speeds 
when the retreating blade is unable to produce the lift 
required to maintain equilibrium with the advancing blade. 
This phenomena is due to the great difference in local 
velocity over the airfoil. The advancing blade encounters 
the velocity of V tip plus the forward velocity. The 
retreating blade encounters V tjp minus forward velocity. At 
high forward speeds, the velocity over the retreating blade 
becomes very small, requiring a very large angle of attack 
to produce the needed lift for equilibrium. Also, the flow 
over the inboard part of the retreating blade is reversed, 
flowing from trailing edge to leading edge. As speed of the 
helicopter increases this reversed flow region moves out 
radially along the blade. Typically at cruise speed as much 
as 30 percent of the blade is experiencing reversed flow. 

The total angle of flow can be determined through either 
of two equations. The first equation is based on the 
average lift coefficient, which is lift curve slope (a) 
times the angle of attack, which will be expressed in terms 
of blade loading (C T /oj . The second equation is based on 
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the flapping angle of the blade. Both of these equations 
are functions of collective pitch angle, 0 Q , longitudinal 
cyclic pitch angle, 0 2 , geometric angle of twist, 0 T , and 
inflow ratio, X. For each of these two equations, the four 
parameters are multiplied by coefficients which are 
functions of the tip loss factor (B) , and the advance ratio 
(;u) . The relationship for the first equation is given by 
Eqn. 25 with the accompanying coefficients given in Eqn. 26, 
27, 28 and 29. 



= xr 1+ 0 o r 2+ 0 r T 3+ 0 2 r 4 


(25) 


T x = 0.5 (B 2 + 0.5^ 2 ) 


(26) 


T 2 = 0.33 {B 3 +0.5\l 2 B) 


(27) 


T 3 = 0.25 ( B 2 (B 2 + \l 2 ) ) 


(28) 


T 4 = 0 . 5|i (S 2 +0 . 25|i. 2 ) 


(29) 



The second relationship is for the longitudinal flapping 
coefficient, a 1 , which is assumed to be zero or very close 
to it. The second set of relationships are given in Eqn. 

30, 31, 32, 33, and 34. 

a-i - 0 = XA 11 +0 o A 12 + 0 r A 13 + 0 2 ^ 14 (30) 
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( 31 ) 



ni 



B 2 (B 2 - 0 . 5\i 2 ) 



12 



8tig 



3 (B 2 - 0 . 5(i 2 ) 



( 32 ) 



A,, = 



2| xB^ 



B 2 - 0 . 5[i 2 



( 33 ) 



*14 



= B 2 + l . 5|i ; 
B 2 - 0 . 5|i 2 



( 34 ) 



Solving the two equations simultaneously in terms of the 
eight coefficients, the collective pitch angle and the 
longitudinal cyclic pitch angle can be determined. The 
resulting equations are as follows: 



0o 



07^13 + 02-^14 + ^11 

Al2 



( 35 ) 



+6 t (A 13 T 2 -A 12 T 3 ) (3g) 

Ai*TrA^T 2 

Using these equations, approximations for the angle of 
attack of the advancing and retreating blades can be 
calculated. The equations for the angle of attack are given 
by Eqn . 37. 
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(37) 



a 



27 0 



0 o- 0 2 +0 r + 



X 

1+Ji 



a 



90 



= 0 o +0 2 +0 r + 



X 

l+n 



Appendix E contains the spreadsheet that was used to 
calculate the angle of attack for three different blade 
twist conditions, 0, -4, and -8 degrees of twist. The minus 
sign indicates that the twist angle decreases from the tip 
of the blade to the root. Fig. 33, 34, and 35 show the 
results of these calculations. 
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Figure 33 Retreating Blade Stall, 0 Twist 
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Figure 34 



Retreating Blade Stall, -4 deg Twist 
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Figure 35 Retreating Blade Stall, -8 deg Twist 
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According to Ref. 8 the blade reaches its operational limit 
when the angle of attack of the tip of the retreating blade 
exceeds the stall angle of attack (AOA) of the airfoil 
section by four degrees [Ref. 8:p. 258]. For the NACA 0012, 
this AOA would be 16 degrees. There are three areas of 
interest on these figures: airspeeds which require less than 
12 degrees AOA where no stall will be encountered, airspeeds 
which require between 12 degrees and 16 degrees AOA where a 
moderate amount of stall is present, and airspeeds which 
require 16 degrees or greater AOA where the severity of the 
stall prohibits flight. Two lines are annotated on the 
above graphs to show the velocity for this stall area. The 
limiting velocities for the three conditions are 64 & 74, 66 
& 76, and 68 & 77 knots respectively. The result agrees 
with theory, because the increase in twist should increase 
the stall envelope, which the figures demonstrate. 
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VIII. 



MODIFICATIONS 



A. LANDING GEAR 

The rear landing gear of the Hummingbird has a narrow 
tread width which is considered to be too narrow in view of 
takeoff and landing stability. It also has no shock 
absorption capability, which is critical for avoiding ground 
resonance. A study was conducted to determine the most 
advantageous route to take to overcome this challenge. The 
options were to leave it as it is, replace it with a skid 
configuration, or a strut configuration. The landing gear 
was redesigned with a strut configuration to overcome both 
of these problems. The tread width was widened from 17 
inches to 27 inches and an oleo was incorporated. The 
spring constant of the oleo was chosen so that the 
Hummingbird could survive a fall from ten feet. Ten feet is 
approximately one rotor diameter, the height to hover out of 
ground effect. 

B. DRIVE TRAIN MODIFICATIONS 

1. Engine to Main Rotor Gear Ratio 

Because of the requirement for the main rotor blades 
to obtain a prescribed Reynolds number range, the rotor rpm 
had to be increased to achieve the tip velocity required. 
Before the gearing ratio was chosen, the decision to derate 
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the engine had to be made. It was shown earlier that for 
the Hummingbird's range of operation, there would always be 
excess power available. Based on this fact and the desire 
to extend the life of the engine, it is recommended that the 
engine be rated to operate at 18 HP vice 25 HP. Fig. 32 
shows that for this horsepower an engine rpm of 5000 is 
required. The present gear ratio is 10:1 and would have to 
be lowered to 6.6:1 to achieve a rotor rpm of 759. 

2. Freewheeling Unit 

The Hummingbird, as stated earlier, has no ability 
to autorotate. Because of this operating limitation, the 
Hummingbird is at great risk of loss to an inflight engine 
failure. It was decided that the most effective way to 
solve this safety hazard was to incorporate a one-way- 

bearing in the gear that would be changed for the gear ratio 

change. After examination of the drive train, the last gear 
in the drive train was selected as the most effective point 
to make the modification. This point was selected for two 
reasons. First, because it was the last gear in the drive 
train, the rotor would have the least mechanical resistance 
during an autorotation. This point is critical because 
every rpm gives that much more cushion in the touch down 

phase of an autorotation, thus increasing the chances of a 

successful autorotation. Second, this gear has easy 
accessibility and would simplify the modification process. 
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IX. 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The goal of this research effort was to design a 
remotely piloted helicopter (RPH) with the capabilities to 
meet the expanding needs of the Department of Aeronautics 
and Astronautics helicopter flight test program. These 
needs have grown rapidly over the past two years. J. G. 
Scott stated in his thesis that the "requirements were: (1) 

a four-bladed rotor head; (2) payload capacity of 
approximately 15 pounds; and (3) a total system cost of no 
greater than $10,000." [Ref. 12] Numerous research 
requirements have arisen since that time, which have 
expanded the envelope of RPH requirements. Lt. M. Borno's 
quarter scale NOTAR research requires a platform of greater 
size and power than is available with any of the 
department's current assets. Continuing HHC research 
requires an RPH that will produce blade Reynolds numbers in 
the range of 1.5 to 1.8 million in order to validate HHC 
performance enhancements. Based on these two research 
requirements and the hindsight of how quickly these 
requirements change, it was imperative to develop an RPH 
that could meet the existing needs and would be flexible 
enough to meet future needs . 
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Upon completion of the initial design phase it was 
determined that two paths were available: (1) complete the 

detailed design in house; or (2) search the open market to 
find an RPH capable of meeting the existing requirements. 
The second path was chosen and resulted in procurement of 
the Hummingbird. The Hummingbird possesses the capabilities 
to meet existing and future requirements with some 
modifications. The NOTAR requirements are fulfilled in the 
Hummingbird's present state, but the HHC requirement 
dictates some modification. These modifications include a 
new three blade main rotor hub to be used with the existing 
rotor blades and a change in the engine to main rotor 
gearing ratio. Two other modifications that would be of 
great value from a safety standpoint are incorporating a 
freewheeling unit in the engine drive system and replacing 
the existing landing gear. 

The performance analysis shows that with the main rotor 
hub change, the Hummingbird will be able to achieve a 
Reynolds number within the range required for aerodynamic 
performance comparison with the 0H-6A. This change will 
reduce the Figure of Merit, but due to low disk loading, 
will have little effect on the Hummingbirds power loading. 

The analysis of the Hummingbird also shows that it will 
have an excess power of nearly 10 HP for HOGE and therefore 
can be operated at a derated power. This would reduce wear 
on the engine and thus enhance its life. The decision to 
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derate the engine should be made prior to changing the 
engine to main rotor gearing, as it will effect the 
operating rpm of the engine. 

The design improvements; modified landing gear and a 
freewheeling unit, will provide the capability to safely and 
successfully recover the Hummingbird from nearly any 
malfunction over a large flight regime. 

The overall analysis of the Hummingbird is that it is an 
ideal platform for subsequent research. It has the 
capabilities to meet the existing needs and the flexibility 
to meet any foreseeable future needs. 

B. RECOMMENDATIONS 

There is almost unlimited potential for the 
Hummingbird's use in further RPH research. Two items that 
need to be undertaken as soon a possible are the design and 
manufacture of a three blade hub and incorporating a 
freewheeling unit in the engine drive train. The first is 
necessary before any HHC work can be undertaken on the 
Hummingbird. Measurement could then be started on the 
control system freeplay and torsional constants. The ground 
work is laid out in Ref. 12 and 13. The second item is 
primarily for safety reasons and to help protect the 
investment of the department. 

The design of a new main rotor hub could be a thesis in 
itself. A great source of expertise that could be of help 
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is Mr. Art Phelps of the Aerostructures Directorate at NASA 
Langley. His directorate owns an RPH identical to the 
Hummingbird, and it would be of value to develop the bonds 
of cooperation that have been established with them. 

Two other areas where valuable research could be done is 
in the design of an improved main rotor blade and the 
incorporation of a new landing gear design. An improved 
main rotor blade should be designed to optimize airfoil 
cross section and blade twist. The manufacture of new 
landing gear should be a moderately easy job which would 
greatly enhance the survivability of the Hummingbird. 

The final recommendation is made in light of all the 
current and possible future research that would hinge on 
the Hummingbird. It is highly recommended that a second 
Hummingbird be purchased. The reason for this is that 
concurrent research on the Hummingbird will inevitably 
require alteration that would adversely affect one or both 
of the areas of research. A case in point is the hub 
redesign and how it will affect the NOTAR research. The 
NOTAR research demands that the rotor operating parameters 
be constant. This demand would not be feasible in the 
tradeoffs that occur in the design process. Therefore, 
having a second helicopter would allow for numerous research 
projects to be accomplished concurrently without adversely 
affecting one another. 
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APPENDIX As MYKLESTAD FORTRAN CODE 



C THIS PROGRAM CALCULATES MYKLESTAD DETERMINANT FOR 
C DETERMINING THE NATURAL FREQUENCY OF A ROTOR BLADE 
C 

C THIS PROGRAM IS FOR AN ARTICULATED HUB ONLY 
C 

C THE PROGRAM IS SET UP TO TAKE A BLADE OF 17 SEGMENTS. 



IMPLICIT DOUBLE PRECISION ( A-H , L, M , O-Z) 

DOUBLE PRECISION T ( 1 : 1 7 ) , SHEAR ( 1 : 17) , MOM ( 1 : 17) , SLOPE ( 1 : 17 ) , 
% Y ( 1 : 1 7 ) , R ( 1 : 1 7 ) , WT ( 1 : 1 7 ) , DELTA1 , XI (1 : 17) ( NR , NRSQ , F , WSQ , 

% L ( 1 : 17 ) ,W,S1,S2,S3,Y1,Y2,Y3,Y4, MOM1 , MOM2 , BS , BM , DELTA 2 ( 
*BSLOPE,BY,AS, AM, ASLOPE , AY, XPLOT ( 1 : 17) , YPLOT(l: 17) 

Q************************************************************* 

C VARIABLE DEFINITION 
C 

C T - CENTRIFUGAL FORCE 

C SHEAR - SHEAR 

C MOM - MOMENT 

C SLOPE - SLOPE 

C Y - DEFLECTION 

C* ************************************************************ 



THIS OPENS THE DATA FILE FOR OUTPUT "WOODPTS . DAT” 

OPEN ( UNIT=3 5 , FILE= ' WOODPTS . DAT' , STATUS^ ' UNKNOWN ' ) 

C 

C************************************************************** 

c 

C THESE ARE THE INITIAL PARAMETERS TAKEN FROM THE BLADE DATA 
C THEY ARE FROM THE BLADE TIP TO THE BLADE ROOR 
C 

C SEGMENT SPACING (IN) 

DATA L/6. 25, 16*3.0/ 

C SEGMENT RADIUS (IN) 

DATA R/ 59. 1,54. 5, 51. 5, 48. 5, 45. 5, 42. 5, 39. 5, 36. 5, 33. 8, 30. 5, 

% 27. 5, 24 .5,21.5, 18.5, 15.5, 12.5,9.5/ 

C SEGMENT WEIGHT (LBS) 

DATA WT/. 459, 13 *.22, .708, .656, .742/ 

C SEGMENT FLAPWISE MOMENT OF INERTIA, Ixx (IN A 4) 

DATA XI/14*.097, .45, .417, .808/ 

Q*************************************************************** 

C ALLOWS FOR DIFFERENT NR INPUTS WITHOUT RECOMPILING 
C 

PRINT *, ' INPUT ROTOR RPM OMEGA ' 

READ *, NR 
NRSQ=NR*NR 

E= 1305000 I MODULUS OF ELASTICITY OF SPRUCE WOOD 

C * *************************************************************** 

C THIS LOOP FOR F MAY NEED TO BE MODIFIED FOR THE LOWER FREQUENCY 
C RUNS AS YOU DO NOT OBTAIN ENOUGH DATA 
C 

DO 400 F= . 1, 30, . 1 
W= ( NR* F ) 

WSQ=W*W 

1 = 1.0 
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c***************************************************************** 

c 

C INITIAL CONDITIONS 
C 

SLOPE ( 1 ) =1 
Y(1)=0 

020 SHEAR ( 1 ) =0 
MOM ( 1 ) =0 
C 

C****************************************************************** 

c 

C MYKLESTAD METHOD 

T ( 1 ) = WT ( 1 ) * N RSQ * R ( 1 ) / (32. 174*12) 

DO 30 J=2 , 17 

T ( J) =T ( J-l ) + ( ( WT ( J) *NRSQ*R ( J) )/(32.174*12)) 

Sl=SLOPE(J-l) * (l+( (T(J-l) *L(J) *L(J) ) / (2*E*XI (J) ) ) ) 

S2=MOM(J-l) *L ( J) / ( E*XI ( J) ) 

S3=SHEAR ( J-l ) *L (J) *L(J) / ( 2 *E*XI ( J) ) 

SLOPE ( J) =S1-S2-S3 
Y l=SLOPE ( J ) *L ( J) 

Y2=T ( J- 1 ) * SLOPE ( J- 1 ) * (L(J) **3) / (3*E*XI (J) ) 

Y 3=MOM ( J- 1 ) * (L(J) **2) / (2*E*XI (J) ) 

Y4 =SHEAR (J-l ) * (L(J) **3) / (3*E*XI (J) ) 

Y ( J) = Y ( J- 1) -Y1+Y2-Y3-Y4 

SHEAR ( J) =SHEAR ( J- 1 ) + ( ( WT ( J ) *WSQ* Y ( J) )/( 32. 174*12)) 

MOM 1=S HEAR ( J- 1 ) *L(J) 

MOM2=T ( J-l ) *(Y(J-1) -Y(J) ) 

MOM ( J ) =MOM ( J - 1 ) +MOM 1 - MOM 2 
030 CONTINUE 

I F ( I . EQ . 2 ) GOTO 200 
BS=SH EAR (17) 

BM=MOM( 17) 

BSLOPE=SLOPE( 17) 

BY=Y (17) 

1 = 2 
C 

C***************************************************************** 

C NEW INITIAL CONDITIONS FOR SECOND COEFFICIENTS 
S LO P E ( 1 ) = 0 
Y(l)=l 
C 

Cl***************************************************************** 

GOTO 020 

200 AS=SHEAR (17) 

AM=MOM (17) 

ASLOPE=SLOPE ( 17 ) 

AY=Y ( 17) 

C 

C DELTA1 VARIABLE IS FOR THE ARTICULATED ROTOR 
C DELTA2 IS FOR THE HINGELESS ROTOR 
C 

C WHERE DELTA GOES TO ZERO IS THE NATURAL FREQUENCY 
C 

DELTA 1= (AM* BY) - (BM*AY) 

DELTA 2= (ASLOPE*BY) - ( BSLOPE* AY ) 

C 

50 FORMAT (IX, 1PE13 . 6 , IX, 1PE12 . 5 f IX, 1PE12 . 5, IX, 1PE12 . 5 , IX , 1 PE 12 . 5 , 
% IX , 1 PEI 2 . 5 ) 
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51 FORMAT (1X,F10.3,2X,F10.3) 

WRITE ( * , 50) W , DELTA 1 , AM, AY , BM , BY 
WRITE (35,50) W , DELTA1 , AM, AY , BM , BY 
400 CONTINUE 
C 

STOP 

END 
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APPENDIX B 



MYKLESTAD NATURAL FREQUENCIES 



First Natural Frequency 
Operating rpm 759 (79.4118 rad/s) 



Delta 



Am 



Coefficients 
Ay Bm 



Frequency 

(rad/s) 



1 . 07 1 52 2 E + 02 
1 . 07 1 52 2 E + 02 

1 . 07 1 522 E + 02 
1 .071522E+02 

1 . 07 1 522 E+02 
1 . 07 1 522E+02 

1 . 07 1 52 2E + 02 
1.07 1522 E+02 
1.071 52 2E+02 
1. 071522E+02 

1 . 07 1 52 2E + 02 

1 . 07 1522 E + 02 



107 . 1522 
(1022 rpm) 



6.82 39 2 E-0 1 
5. 74431E-01 
4 . 66470E-01 
3.58510E-01 
2. 50549E-01 
1.42588E-01 
3 . 46273E-02 
-7.33 3 35E-02 
-1.8 1294E-01 
-2 . 89255E-01 
-3 .972 16E-01 
-5.0517 7 E-0 1 



0.0 



1. 79121 E+04 
1 . 79 1 2 1 E+04 
1.79121 E+04 
1 . 79121E+04 
1 . 79 12 1E+04 
1.7912 1E+04 
1 . 79121E+04 
1 . 79 122E+04 
1 . 79122E+04 
1 . 79122E+04 
1 . 79122E+04 
1.79122E+04 



= 0 

1.7912 1 E+04 



7 . 48559E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 4 8560E + 00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48560E+00 
7 . 48 561E + 00 



7.48560 



-1.12 3 8 4 E+ 06 
-1 . 12384E+06 
-1 . 12384 E + 06 
-1 . 12384E+06 
-1. 12384E+06 
-1. 12384E+06 
-1. 12384E+06 
-1. 12 384 E + 06 
-1. 12 3 84 E + 06 
-1. 12384E+06 
-1 . 12384E+06 
-1. 12384E+06 



-1 . 12 384 E + 06 



Frequency at which Delta 



By 



-4 . 69659 E+02 
-4 . 69659E+02 
- 4 . 69 659 E+02 
-4 . 69 659E + 02 
-4 . 69 659E+02 
-4 . 69659E+02 
-4 . 69659E+02 
-4 . 69659 E+02 
-4 . 69659E+02 
-4 . 69659E+02 
-4 . 69659E+02 
-4 . 69659E+02 



-4 . 69659E+02 
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2nd Natural Frequency 
Operating rpra 759 (79.4118 rad/s) 



Frequency 

(rad/s) 



3 . 38 1635E+02 

3 . 38163 5E + 02 
3.38 1636E+02 

3 . 38 163 7E + 02 
3 . 381638E+02 
3 . 381639E+02 
3.38 1639E+02 

3.38 164 0E + 02 
3 . 38164 1E+02 
3.331642Ef02 
3 . 3 8 1 64 3 E + 02 

3.38164 3 E + 02 

3.38164 4 E + 02 

3 . 38164 5E + 02 

3.38164 6E + 02 
3 . 38 1647E+02 

3.38 164 7E + 02 
3. 381648Ef02 

3.38 1 64 9 E + 02 



Frequency at 
338 . 1641 
(3229 rpm) 



Coefficients 



Am 



Ay 



Bm 



By 



Delta 



-2 . 158 18E+01 
-1 . 90265E+01 
-1 . 6471 1E+01 
-1 . 3 9158 E + 01 
-1 . 13605E+01 
-8 . 80514E+00 
-6 . 24980E+00 
- 3 . 694 4 6E + 00 
-1. 13912E+00 
1.41623 E + 00 
3 . 97157 E + 00 
6 . 5 2 692 E + 00 
9.08227E+00 
1 . 16376E+01 
1. 41930E+01 
1 . 67483 E + 01 
1.93037E+01 
2. 18591E+01 
2.4414 4 E + 0 1 



whci Delta 

0.0 



3 . 50008 E + 05 

3 . 50008 E + 05 

3 . 50008 E + 05 
3 . 50009E+05 
3 . 50009E+05 
3 . 50009E+05 
3 . 50010E+05 
3 . 50010E+05 
3 . 50010E+05 
3 . 500 10E+05 
3 . 5001 1E+05 
3 . 5001 1 E+05 
3.5001 1E+05 
3 . 5001 1E+05 
3 . 5001 2E+05 
3 . 5001 2E+05 
3 . 50012E+03 
3 . 5001 2E+05 
3 . 50013 E + 05 



3 . 50010E+05 



8 . 05647E+01 

8 . 0564 7E + 01 
8 . 05647Ef 01 
8 . 05648E+01 
8. 05648E+01 

8 . 0564 9 E + 0 1 

8 . 0564 9E + 01 
8 . 05650E+01 
8 . 05650E+01 
8 . 05651E+01 
8. 05651E+01 
8 . 05652E+01 
8 . 05652E+01 
8 . 05652E+01 

8 . 0565 3 E+01 

8 . 05653 E + 01 

8 . 05654 E + 0 1 

8 . 05654 E + 01 
8 . 05655E+01 



8 . 05650E+0 1 



-5. 16530E+06 
-5. 1 65 3 0E + 06 
-5. 16530E+06 
-5. 16531E+06 
-5. 1 65 3 1 E 1-06 
-5. 16531E+06 
-5. 16532E+06 
-5. 16532E+06 
-5. 16532E * 06 
-5. 16532E+06 
-5. 16533E+06 
-5. 16533E+06 
-5. 16533E+06 
-5. 16534E+06 
-5. 1653 4E + 06 
-5. 16534E+06 
-5. 1653 4 E + 06 
-5. 1653 5E+06 
-5. 16535E+06 



-5. 16532E+06 



-1. 18895E+03 
-1. 18895E+03 
-1. 18895E+03 
-1 . 18895Ef03 
-1. 18895E+03 
-1. 18895E+ 03 
-1. 18895E+03 
-1 . 18895E+03 
-1. 18895E+03 
-1. 18895E+03 
-1 . 18895E+03 
-1 . 18895E+03 
-1. 18895E+03 
-1. 18895E+03 
-1. 18895E+03 
-1.18895E+03 
-1. 18895E+03 
-1 . 18895E+03 
-1 . 18895E+03 



-1. 1 8895E+03 



= 0 
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3rd Natural Frequency 
Operating rpm 759 (79.4118 rad/s) 



Frequency 
( rad/s) 



7 . 5 8 3 7 2 5 E+ 02 
7 . 58 3 72 5E + 02 
7.58372 5E + 02 
7.583725Ef02 
7. 583725E+02 
7 . 583725E+02 
7.583725Ef02 
7.583725E+02 
7 . 583725E+02 
7 . 583 725E + 02 
7 . 583725E+02 
7. 583725E+02 
7 . 583726Ef 02 
7. 583726E+02 
7. 583726E+02 
7 . 583726E+02 
7. 583726E+02 
7 . 58 3 72 6E+ 02 



Frequency at 
758.3725 
(7242 rpm) 



Delta 



2 . 20479E+01 
2 . 00033E+01 
1 . 79588E+01 
1 . 59 1 4 3 E + 0 1 
1.38698E+01 
1. 1 8 2 5 3 E + 0 1 
9 . 78 079 E + 00 
7.73628E+00 
5 . 69 177E ^00 
3 . 6 4 7 2 5E + 00 
1. 60272E+00 
-4 . 4 17 68E-01 
-2 . 4 8627E + 00 
-4 . 5308 1 E + 00 
-6 . 575 3 2 E+00 
-8 . 61982Ef 00 
-1 . 06643E+01 
-1. 27089E+01 



which Delta 

0. 0 



Am 



8.24037E+06 
8 . 24038E+06 

8 . 2 4038 E+ 06 
8 . 24038E+06 
8.24038Et06 
8. 24038E+06 

8 . 2 4 03 8 E + 06 

8 . 2 4 03 8E+ 06 

8 . 2 4 038 E+06 

8 . 2 4 03 8 E + 06 
8 . 24038E+06 
8 . 24038Et06 
8 . 24038E+06 

8 . 2 4 038 E+06 
8. 24038E+06 
8 . 24038Ef06 
8 . 24038E+06 
8 . 24038E+06 



0 

8.24038Ef06 



Coe f f i 
Ay 



8 . 12 28 5E + 02 
8. 12285E+02 
8 . 12286E+02 

8 . 12 286E + 02 
8 . 12 286E + 02 
8. 12286E+02 
8 . 122 86E + 02 
8 * 1 2286E+02 
8. 12286E+02 
8. 12286E+02 
8. 122 86E + 02 
8. 12286E+02 
8. 12286E+02 
8 . 12286E+02 
8.12286E+02 
8. 12286Ef02 
8. 12286E+02 

8 . 1 2 2 8 6 E + 02 



8 . 12286E+02 



ients 

Bm 



-7.04041Ef07 
-7 . 0404 1E + 07 
-7.04041E+07 
-7.04041E+07 
-7.04041E+07 
-7 . 04 04 1 E + 07 
-7.04041E+07 
-7 . 04041E+07 
-7 . 04041E+07 
-7.0404 1 E + 07 
-7 . 04041E+07 
-7 . 0404 1E + 07 
-7 . 04 04 1 E + 07 
-7 . 04 04 1 E + 07 
-7 . 04 04 1 E + 07 
-7 . 04 04 1 E + 07 
-7. 04041 E+07 
-7.04041E+07 



-7 . 04 04 1 E + 07 



By 



-6. 94 OOOE + 03 
-6 . 94 000E*03 
-6.94 OOOE + 03 
-6 . 94 000E + 03 
-6 . 94 OOOE+03 
-6.9 4 000E + 03 
-6.9 4 000E+03 
-6.94000E+03 
-6 . 9 4 000E+03 
-6 . 94000E+03 
-6. 94000Ef03 
-6.9 4 000E+03 
-6 . 94 000E+03 
-6 . 9 4 00OE + 03 
-6 . 94 000E + 03 
-6.94000Ef03 
-6 . 94 OOOE + 03 
-6.94 OOOE + 03 



-6.94 000E + 03 



82 



4th Natural Frequency 
Operating rpra 759 (79.4118 rad/s) 



Frequency 

(rad/s) 



1. 448090E + 03 
1. 448091E+03 
1.448092E+03 
1. 448092E + 03 
1. 448093E + 03 
1. 448094E+03 
1. 448095E+03 
1 . 4 4 8096E + 0 3 
1. 448096E + 03 
1. 448097E + 03 
1. 448098E + 03 
1. 448099E+03 
1. 448100E + 03 
1. 448100E+03 
1.448 101E+03 
1.448 102E+03 
1. 448103E + 03 
1. 4 4 8 104 E + 0 3 



Delta 



-4 . 21455E+04 
-3. 88022E+04 
-3 . 54 589E+04 
-3 . 21156E+04 
-2 . 87722E+04 
-2 . 54 288E + 04 
-2 . 208 54 E+04 
-1.87 42 0E+04 
-1 . 53985E+04 
-1 . 2 0551E + 04 
-8 . 71160E+03 
-5 . 36809E+03 
-2. 02456E+03 
1 . 31898E+03 
4 . 66255E+03 
8 . 006 13E+03 
1. 13497E+04 
1. 46934E+04 



Am 



2 . 7 7 806 E + 08 
2 . 77807E+08 
2 . 77808 E+08 
2.77 8 09 E + 08 
2.778 10E+08 
2 . 778 11E+08 
2 . 778 12E+08 
2 . 778 13E+08 
2 . 778 14 E + 08 
2 . 778 15E+08 
2 . 778 16E+08 
2 . 778 17E+08 
2 . 778 18E+08 
2 . 77$ 19E+08 
2 . 77820E+08 
2 . 7782 1E+08 
2 . 77 822E + 08 
2 . 7 7 8 2 2 E+ 08 



Frequency at which Delta = 0 



1448.100 0.0 2 . 77 8 19 E + 08 



Coefficients 

Ay Bra By 



1 . 184 96E + 04 -1 
1 . 18497 E+04 -1 
1.18497E+04 -1 
1 . 18497 E + 04 -1 
1 . 18498E+04 -1 
1.18498E+04 -1 
1 . 18498E+04 -1 
1.18499E+04 -1 
1.18499E+04 -1 
1. 184 99E + 04 -1 
1 . 18499E+04 -1 
1 . 18 500E+ 04 -1 
1 . 18500E+04 -1 
1 . 18500E+04 -1 
1 . 1850 1E+04 -1 
1.18501E+04 -1 
1 . 18 50 1E + 04 -1 
1 . 18 502 E+04 -1 



. 7 52 50E ♦ 09 -7 
.75251E+09 -7 
.75251E+09 -7 
. 75252E+09 -7 
. 7 5253E+09 -7 
.75253E+09 -7 
. 7 5254 E +09 -7 
. 7 52 54 E + 09 -7 
. 75255E+09 -7 
. 7 52 55E +09 -7 
. 75256E+09 -7 
. 7 5257 E + 09 -7 
. 75257E+09 -7 
. 75258E+09 -7 
. 75258E+09 -7 
. 75259E+09 -7 
. 7 52 59E + 09 -7 
. 7 52 60E+09 -7 



1 . 18500E+04 -1.752 58 E + 09 -7 



(13828 rpm) 



4 7 5 18 E+04 
4 7520E + 04 
4 752 1E + 04 
4 752 3E + 04 
4 7 52 5E + 04 
47527E+04 
4 7529E+04 
47531E+04 
4 7 532 E + 04 
47534 E* 04 
4 7536E + 04 
47538E+04 
4 754 0E+04 
47542E+04 
4 7 544 E+04 
4 754 5E+04 
4 754 7 E + 04 
4 754 9E+04 



4754 1E+04 
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First Natural Frequency (Light Blade) 



Operating rpm 759 (79.4118 rad/s) 



Frequency Delta 

(rad/s) 

Am 



Coefficients 
Ay Bm 



9 . 987550E+0 1 
9.987551E+01 
9. 987552E+01 
9 . 98 7553E + 0 1 
9 . 98 7 553 E+0 1 
9 . 98 7 554 E+ 0 1 
9. 987 555E + 0 1 
9. 987556E+01 
9 . 98 7 557 E+0 1 
9 . 987557E+01 
9 . 987558 E + 0 1 
9 . 987559E+0 1 
9 . 987560E+01 



3 . 7 5536E-01 
3 . 32246E-01 
2 . 88956E-01 
2 . 45666E-01 
2 . 02376 E-01 
1.59087E-01 
1 . 15797E-01 
7.25066E-02 
2.92 1 66E-02 
-1 . 4 07 3 3E-02 
-5. 73633E-02 
-1 . 00653E-0 1 
-1 . 43943E-01 



6.90470E103 
6 . 90470E+03 
6 . 904 7 0E + 03 
6 . 904 70E+03 
6 . 90470E+03 
6 . 9047 1E+03 
6.9047 1E+03 
6 . 9047 1E+03 
6 . 9 04 7 1 E + 0 3 
6.90471E+03 
6.9047 1E+03 
6.90471E+03 
6.90471Et03 



2 . 51694E+00 -4.13426E+05 
2 . 5 1 69 4 E+ 00 -4.13426E + 05 
2 . 51694 E+00 -4.13426E+05 
2.51694 E+00 -4.13426E + 05 
2 . 51694 E+00 -4.13426E+05 
2 . 51694E+00 -4.13426E+05 
2 . 51694E+00 -4.13426E+05 
2 . 51694E+00 -4.13426E+05 
2 . 51694 E + 00 -4.13426E + 05 
2 . 5 1695E+00 -4.13426E+05 
2 . 5 169 5E + 00 -4.13426E + 05 
2 . 5 1695E+00 -4.13426E+05 
2 . 5 1695E+00 -4.13426E+05 



Frequency at which Delta = 0 

99.87557 0.0 6.90471E+03 2.51695 -4.13426E+05 



(954 rpm) 



By 



- 1 . 507 04 El02 
-1 . 50704 E + 02 
-1. 50704E+02 
- 1 . 50704 E+02 
-1. 50704E+02 
-1 . 50704E+02 
-1.50704E+02 
- 1 . 507 04 E+02 
- 1 . 50704 E + 02 
- 1 . 507 04 E+02 
-1- 50704E+02 
- 1 . 50704 E+02 
-1. 50704E+02 



-1. 50504E+02 
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2nd Natural Frequency (Light Blade) 
Operating rpra 759 (79.4118 rad/s) 



Frequency 

(rad/s) 



4 . 004 9 6 1 E+ 02 
4.004961E+02 
4 . 004961E+02 

4 . 004 9 6 1 E+02 
4 . 004961E+02 
4.G04962E+02 
4 . 004962E+02 

4 . 004 9 62 E+ 02 
4 . 004962E+02 
4 . 004962E+02 

4 . 004 962E + 02 
4 . 0049 62E ^02 

4 . 004 9 62 E + 02 
4 . 0049 62E + 02 

4 . 004 9 6 2 E + 02 
4 . 004962E+02 
4. 004962 E+02 
4 . 004962E+02 



Frequency at 
400.4962 
(3824 rpm) 



Delta 



-1 . 03982E+00 
-8.98344 E-0 1 
-7 .56864 E-01 
-6 . 15385 E-01 
-4.73 906E-01 
-3 . 3 2 4 2 6 E-01 
-1.90947E-01 
-4 .9 4 674 E-02 
9 . 20121E-02 
2.3349 2 E-0 1 
3 . 74971E-01 
5 . 16451E-01 
6.579 3 OE-0 1 
7.994 10E-01 

9 . 4 0889E-01 
1 . 08237E+00 
1.2238 5 E + 00 
1 . 36533E+00 



which Delta 

0.0 



Am 



1.91468E+05 
1 . 91468 E+05 
1.914 68E + 05 
1 .91468E+05 
1 .91468E+05 
1.91468E+05 
1 . 91468E+05 
1 . 91468E+05 
1 .91468E+05 
1 . 9 1468E+05 
1 . 9 1 4 68E + 05 
1 .91468Ef05 
1 . 9 1468E+05 
1 .91468E+05 
1 .91468E+05 
1 .91468E+05 
1.91468E+05 
1 .91468E+05 



= 0 

1 . 9 1 4 68 E+05 



Coef t ici 
Ay 



3 . 053 55E + 01 

3 . 053 55E + 0 1 

3 . 053 55E+01 
3 . 05355E+01 

3 . 053 55E+01 

3.053 5 5E + 0 1 

3 . 053 55E + 01 
3 . 05355E+01 
3.05356E+01 
3 . 05356E+01 

3 . 053 56E + 01 
3 . 05356E+01 
3 . 05356E+01 

3 . 053 56E+01 

3 . 053 56E + 01 

3 . 053 56E + 01 
3 . 05356E+01 

3 . 053 56E + 01 



3 . 053 56E + 01 



Bra 



-2.9198 1E+06 
-2.91981 E+06 
-2.9198 1 E+06 
-2.91981 E+06 
-2.91981E+06 
-2.91981E+06 
-2.91981E+06 
-2.9198 1E+06 
-2.91981E+06 
-2.9198 1 E+06 
-2.9198 1E+06 
-2.9198 1 E+06 
-2.9198 1 E+06 
-2 . 9 198 1E+06 
-2. 91982 E+06 
-2 . 9 1 982 E + 06 
-2.9198 2 E + 06 
-2 . 91982E+06 



-2 .9198 1 E+06 



By 



-4 . 65656E+02 
-4 . 65656E ^02 
-4 . 65656E+02 
-4 . 65656E+02 
-4.65656E+02 
-4 . 65656E f02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65 6 56E+02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65656E+02 
-4 . 65656E + 02 
-4.65656E+02 
-4 . 65656E+02 



-4 . 65656E+02 
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3rd Natural Frequency (Light Blade) 
Operating rpm 759 (79.4118 rad/s) 



Frequency Delta 
( rad/s) 



Am 



Coef t icients 
Ay Bm 



By 



1 . 062 3 60E + 0 3 
1 . 06236 1E+03 
1 . 06 2 3 62 E t 03 
1.062362E+03 
1.062363 E + 0 3 
1. 062364E+03 
1.06236 5 E + 0 3 
1 . 0623 65E + 03 
1. 062366E+03 
1 . 062367E+03 
1 . 06 2 3 68 E + 03 
1 . 06 2 3 6 9 E + 0 3 
1 . 062369E+03 
1 . 062370E+03 
1.06237 1 E + 03 



Frequency at 
1062 . 366 
(10145 rpm) 



1 . 49707E+03 
1 . 3 1505E+0 3 
1. 13303E+03 
9.510 10E+02 
7 . 68988E+02 
5 . 86965E+02 
4.0494 1E + 02 
2 . 22916E+02 
4 . 08897E+01 
“1.411 37E+02 
-3.23166Ef02 
-5.05195E+02 
-6 . 87 2 2 5E + 02 
-8 . 692 56E+02 
-1.05 129E+03 



which Delta 

0.0 



7.0618 5E + 06 
7. 06187 E+ 06 
7.06189Ef06 
7 . 06 192E+06 
7.06194 E + 06 
7 . 06197E+06 
7 . 06199E + 06 
7.06201 E+06 
7 . 062 04 E+06 
7 . 06206 E+06 
7 . 06209 E+06 
7.06211 E+06 
7 . 062 13E+06 
7 . 062 16E+06 
7.0621 8 E + 06 



= 0 

7 . 06204E+06 



4 . 97993E+02 
4 . 97995E+02 
4 . 97996E+02 
4 . 97997E+02 
4 . 97999E+02 
4 . 98000E+02 
4 . 98001E+02 
4 . 98003E+02 
4 . 98004E+02 
4 . 98005E+02 
4 . 98007E+02 
4 . 98008E +02 
4.98009E+02 
4 .98011 E+02 
4 .98012 E+02 



4 . 98004 E f02 



-6 . 334 27E+07 
-6 . 33429E+07 
-6. 33431E+07 
-6.334 33E + 07 
-6 . 33435E+07 
-6 . 33437E+07 
-6. 33439E+07 
-6.3344 1E+07 
-6 .3344 3E + 07 
-6. 33444E+07 
-6. 33446E+07 
-6 . 3 34 4 8E + 07 
-6. 3 34 50E + 07 
-6 . 33452E+07 
-6. 33454E+07 



-6.33443 E+07 



-4 . 46686E*03 
-4 .46687E + 03 
-4 . 46688EV03 
-4 . 46689E+03 
-4 . 46690E+03 
-4 . 4 6 69 1E+03 
-4 . 4 669 2E+03 
-4 . 46693E+03 
-4 . 46694E+03 
-4.46695Ef03 
-4 . 46696E+03 
-4 . 4 6697E+03 
-4.4 6698E+03 
-4 . 4 6699E+03 
-4.4 6700E+03 



-4 .46694 E + 03 
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4th Natural Frequency (Light Blade) 
Operating rpm 759 (79.4118 rad/s) 



Frequency Delta 

(rad/s) 

Ara 



Coefficients 
Ay Bm 



By 



2 . 0909 1 3 E + 0 3 
2 . 0909 1 3 E *-03 
2.090914E+03 
2 . G9Q915E+03 
2 . 09091 6E+03 
2 . 0909 17E+Q3 
2.090917E+03 
2. 090918Ei^03 
2 . 0909 19E+03 



-1.4764 6E + 04 
-1 . 15527E+04 
-8 . 34 069E + 03 
-5. 12869E+03 
-1 . 91669E+03 
1 . 2953 3E + 03 
4 . 50736E+03 
7.7194 1E+03 
1 . 09 3 1 5E + 04 



2 . 52950E+08 
2.52951E+08 
2.52952E+08 
2 . 52952E+08 
2.52953E+08 
2 . 52953E+08 
2.52954E+08 
2 . 52955E + 08 
2 . 52955E+08 



9. 15755E+03 
9. 15757E+03 
9. 15759E + 03 
9. 15761 E+ 03 
9. 15762E+03 
9 . 1 5764 E + 03 
9. 15766E+03 
9. 15768E+03 
9. 15769E+03 



- 1 . 68 095 E + 09 -6. 
- 1 . 68096E+09 -6. 
- 1 . 68Q96E+G9 -6. 
- 1 . 68Q97E+09 -6. 
-1 . 68097E 109 -6. 
-1 . 68097E+09 -6. 
-1.68098E+09 -6. 
-1 . 68098E+09 -6. 
-1 . 68098E+09 -6. 



Frequency at which Delta = 0 

2090.917 0.0 2 . 5295 3 E + 08 9.15764E+03 -1.68097E+09 -60 



(19967 rpm) 



08 5 55 E+ 04 
08556E+04 
08557 E + 04 
08558Ef04 
08559Ef04 
08560E+04 
08561E+04 
08563 E t04 
08 564 E + 04 



. 8 560E + 04 
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APPENDIX C: MODE SHAPE FORTRAN CODE 



***** 

* This Program is used to calculate the mode shapes once 

* the natural frequencies are determined by the Myklestad 

* determinant method 

* * * * * 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION RSTA ( 2 4 ) , DMASS ( 2 4 ) , El ( 2 4 ) , CENT ( 2 4 ) , SHEAR ( 2 4 ) 
DOUBLE PRECISION DMOM (24) ,SLP(24) , DEFL ( 24 ) , ix ( 2 4 ) , dm ( 2 4 ) , 

%dx ( 2 4 ) , mod, pi , p 

OPEN ( 101 , FI LE= ' MODE . OUT' , STATUS= ' UNKNOWN ' ) 



Pi 


= 


3 . 141592654 


i pi 


g 


= 


32 . 174 


! GRAVITATIONAL CONSTANT 


e 


= 


0. 133 


1 HINGE OFFSET RATIO 


RAD 


= 


68.0 


1 ROTOR RADIUS 


RO 


- 


e*RAD 


1 HINGE POINT 


C 


= 


6.5 


1 CHORD 


to 


= 


.13 


1 ROOT SECTION THICKNESS (t/c) 


tt 


= 


.13 


1 TIP SECTION THICKNESS (t/c) 


P 


= 


. 0 


1 DENSITY OF SPRUCE Ibrn/in'I 


mod 


= 


1305000.0 


1 MODULUS OF ELASTICITY SPRUCE 



PRINT *, 'ENTER ROTOR RPM AND NATURAL FREQUENCY' 

READ *,rpm,wrpm 
RV=rpra*2.0*pi/60.0 
w=wrpm*2 . 0*pi/60 

1 THE FOLLOWING ARRAYS ARE NUMBERED FROM BLADE TIP TO ROOT: 

1 SEGMENT RADIUS (IN) 

DATA RSTA/ 6 5. 6, 62. 0,59. 6, 57. 2, 54. 8, 5 2. 4, 50.0, 47.6,45.2,42.8, 
%40. 4, 38. 0,35. 6, 33. 2, 30. 8, 28. 4, 26. 0,23. 6, 21. 2, 18. 8, 16. 4, 14.0, 
111.6,9.2/ 

I SEGMENT SPACING (IN) 

DATA dx/ 4.8,23*2.4/ 

l SEGMENT WEIGHT (LBS) 

DATA dm/. 3 53, 20*. 177, . 302 , . 252 , . 470/ 

i SEGMENT FLAPWISE MOMENT OF INERTIA, IXX (IN*4) 

DATA ix/21*.478, .48, .400, .845/ 



C DO 10 J=24,l,-1 

C x — (J-l) *dx(LL)+0.5*dx(LL) 

C RSTA ( 2 5-J ) = X 

C t = ( ( RAD-RO-X) * ( tO-tt ) / (RAD-RO) 3 tt) *C 

C dm = (2. 0*6. 60*0. 2850+2. 0*t*0. 1250) *p*dx(LL) 

C DMASS ( 2 5-J ) » dra/g 

C i - (2.0*6.85*. 2850) *( ( .5* (t-. 1425) ) **2.0) 

C El ( 2 5-J ) - i *mod 

CIO CONTINUE 

I CENTRIFUGAL FORCES COMPUTED ALONG THE BLADE 
DMASS ( 1 ) =dm ( 1 ) /g 

CENT ( 1 ) =DMASS ( 1 ) *(R0+RSTA(1) )*(RV**2)/12 
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DO 30 J=2 , 2 4 

DMASS (J ) =dm ( J) /G 
El (J)=ix(J) *mod 

CENT( J)=CENT(J-1) + DMASS (J) * (R0+RSTA(J) ) * (RV**2) /12 
30 CONTINUE 

I ITERATION LOOP 

SLP(l) =0.0 I BC'S AT TIP FOR Ay AND Am 

DEFL(l) = 1.0 
DMO M ( 1 ) =0.0 

SHEAR ( 1 ) = 0.0 
DO 40 KK=1 , 2 

DO 50 I.L=2 , 24 I MYKLESTAD INTEGRATION 

SLP 1=SLP ( LL- 1) * (1+CENT(LL-1) *(dx(LL) * *2) /2/EI (LL) ) 
SLP2=DMOM(LL-l) *dx ( LL) / El ( LL) 

SLP3=SIIEAR (LL-1 ) * (dx(LL) * *2 ) / 2/ El ( LL) 

SLP ( LL) =SLP1-SI.P2-SLP3 
DEFL1=SLP ( LL) *dx(LL) 

DEFL2=CENT ( LL- 1 ) *SLP ( LL-1 ) * (dx ( LL) * *3 ) / 3 / El (LL) 
DEFL3=DMOM ( LL-1 ) * (dx(LL) * *2 ) / 2/EI ( LL) 

DEFL4=SHEAR ( LL- 1 ) * (dx(LL) **3) /3/EI (LL) 

DEFL(LL) =DEFL(LL-1) -DEFL1+DEFL2-DEFL3-DEFL4 
SHEAR (LL)=SHEAR( LL-1) +DMASS ( LL) * (w**2) *DEFL ( LL) / 12 
DM0M1=SHEAR ( LL- 1 ) *dx(LL) 

DMOM2=CENT(LL-l) * (DEFL(LL-l) -DEFL(LL) ) 

DMOM ( LL) =DMOM ( LL- 1 ) +DMOM1-DMOM2 
50 CONTINUE 

IF (KK.EQ.l) THEN i COEF. COMPUTATION 

Ay=DEFL( 24 ) 

Am=DMOM(24) 

ELSE 

By=DEFL ( 24 ) 

Bra=DMOM(24) 

END IF 





SLP(1)-1.0 
DEFL ( 1 ) =0 . 0 


! BC'S 


FOR 


By 


AND 


Bm 


40 


CONTINUE 


1 LOOP 


FOR 


By 


AND 


Bm 



DET=Am*By-Ay*Bm 
PRINT * , DET 



SLP ( 1 ) = -Ay/By 

DEFL(l) = 1.0 
DMOM ( 1 ) = 0.0 

SHEAR (1) =0.0 

DO 60 LL=2,24 ! MYKLESTAD INTEGRATION 

SLP1=SLP ( LL- 1 ) * ( 1+CEHT ( LL-1 ) *(dx(LL)**2)/2/EI(LL) ) 
SLP2=DMOM ( LL-1 ) *dx (LL) / El ( LL) 

SLP3=SHEAR ( LL- 1 ) * (dx(LL) **2) /2/EI (LL) 

SLP ( LL) =SLP1-SLP2-SLP3 
DEFL1=SLP (LL) *dx(LL) 

DEFL2=CENT(LL-1) *SLP (LL-1 ) * (dx (LL) **3) /3/EI (LL) 
DEFL3=DMOM ( LL-1 ) *(dx(LL) **2) /2/EI (LL) 

DEFL4 =SHEAR ( LL-1 ) * (dx(LL) **3) /3/EI (LL) 

DEFL(LL) =DEFL(LL-1) -DEFL1+DEFL2-DEFL3-DEFL4 
SHEAR (LL)=SHEAR( LL-1) 1 DMASS ( LL) * (w**2) *DEFL(LL) / 12 
DMOM 1=SH EAR ( LL-1 ) *dx(LL) 

DMOM2=CENT ( LL-1 ) * (DEFL(LL-l) -DEFL(LL) ) 

DMOM(LL) =DMOM ( LL- 1 ) +DMOM1-DMOM2 
60 CONTINUE 
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DO 100 J=1 , 2 4 

WRITE (101,*) J,RSTA(25-J) ,DEfL(25-J) 
100 CONTINUE 



* 

* This program was written by Lt Matt Avila 

* Modifications were made by Lt J.L. Vandiver 

* 



END 
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APPENDIX D: PRELIMINARY POWER CALCULATIONS 



Thaai© Helicopter Devign 



Preliminary Pow *r Calculation* 

Vatreble© ©nd conatanl© bom Mainrotor worbtheal ink 



8 * 




0 07 


a 


190 00 


Ibe 




t* 




03417 


a vop- 


20* 02 


knole 


430 00 


b- 




3 


vt- 


40 00 


knota 


07 01 tpe 


A* 




100 04 


B~ 3 Afl. 


10 4* 






Itgia* 


000133 


om*g« 


70 41 


ta d/a 


73* rpm 


01- 




1 40 


pa* Ci- 


30OC 03 






81- 




0 03303 


Cl(avg). 


0 30 






Mu- 




0 13003 


rt>©- 


3 30£ 03 










i soe^oa 


(*4000*3 


1 03E 03 






NACAO0I3 


CLmaa Cdo 


Clalph© 


Met 


Alpha (auU) 








1 29 0 0007 


9 73 


0 71 


12 


P©». 


a* reared M Hover ooc 












P!-(T~ 1 3/(2*rho*A) ~ SJ/BASO HP 


0- Tip ioaa lector 












B - 1 (2*Ct) - 9* 












• • 


0 973* 






PI - 


4 093 HP 




BfAJt). 


0 970* Muat correct toe rho In Cl 






PI (ail) - 


9 333 HP 














Uwd new 8 and rho correction 










P© - Cdo*rho*eigma*A*V«p * 3/4400 HP 










Po • 


4 300 HP 








FlgUt 


-e ot Mot* 














FM- ktaai pwt i actual pwr 
FM-P1'(P1 Po) 












FM- 


0 33 


FM range 


7 to 0 












Thla value ia a low due to Mgh tip apeed and 
low Afl *la related to DL end 8L 


Power required to ho ver KJ« (431) 












Pi<1G£)/PI(OG£)-{ 0 !37**(hA3)' 


*4 * 0 7O50*(h/O) “ 3- 1 


43M*(a/Or 2H 34J2*(h/O)40 3147) 








8- height ©4 *>e rotor above Ova eutlace 










a- 3 3 


A 












0- ro»c x diameter 




h/D - 


0 2200 








0- 11 33 


A 












PipG£)/P1(OG£). 


0 7473* 










PI (MIC)- 


3 70 HP 








faraiAt Dfa| ©1 l/U futt/af* 












Equivalent Flat Plata Area 














a- oi 


Thla weed- lav mined by a chert bom a 8rko*M r 8tody 








Paraait© Power (Pp)-O O'rho’to'VI ~ 3 










Pp- 0 93 


HP 




T ©4 rotor power r*q - 
9X*Pi (rotor) at a hover and 








Po(ti*ghl)- Po(hovet)*{1 *4 3*Uu 
P1(fll 9 at). 2/(rt>o , A'Vf) 




3* ft (rotor) In midrange fwd vei 


Vei (Vnola) 


V - CPI 


Panel la Pro Ole 


Induced 




Total Adv Bid Mach 




0 


o oo 


000 4 40 


4 00 


0*4 


10 10 0 40 




20 


33 70 


007 4 91 


293 


094 


7 00 0 43 




40 


07 3t 


0 93 4 03 


t 20 


033 


0*3 04* 




30 


04 30 


1 04 0 00 


1 01 


0 21 


7 33 0 4# 




M 


101 37 


1 00 0 3* 


004 


0 24 


0 34 0 4* 




70 


110 19 


3 00 #70 


0 72 


0 2* 


0 3* 0 31 




00 


139 02 


4 2* #10 


01) 


099 


11 34 0 33 




00 


191 00 


00* 0 3* 


030 


00* 


13 04 0 34 




too 


100 70 


0 33 7 0* 


0 91 


0 7* 


10 00 0 99 


4000 


NH(k«r •> rtvo- 0001017 






Altitude correcOon © 

Adv Blade 






Pereeile 


ProOle bxjuced 


TiWpwr 


Total 


Mach Ho PI (all)- P1(al)*(rho(al)/rho(aJt)) 






000 


3 34 0 10 


0 07 


10 9# 


0 41 Po(aA)- Po(aF)*rho(alt)/rho(d) 






0 09 


3 03 3 13 


04* 


7 30 


0 44 Pp(al1)-Pp{d)»rho(*U)/rho(d) 






0 43 


30* 1 37 


0 2# 


0 l# 


0 47 PI hover weed dlAarent con acton 






004 


40* 1 29 


0 I# 


03* 


0 4$ A Incorporta led change In Cl- > 






1 49 


4 32 1 04 


0 20 


7 02 


0 30 






3 30 


4 00 0 90 


0 23 


003 


033 






343 


4 93 0 70 


0 4* 


0 30 


033 






4 44 


9 20 0 70 


034 


II 41 


035 






0 70 


9 OS 0 93 


0 03 


13 07 


030 
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APPENDIX E 



BLADE TIP PITCH ANGLE WORKSHEET 



Ihaata H#*V htlco^M OtMgA 



Pitch AnqU D+tfmtnotlon. 

ln«ow talk* lam- (| Kp/W W/(»h«>*A*Vf))/VUpl 

Advanca Raao Uu> Vt/Vkp 

P*/a#*la Powa* jT’pl-G 6*rho*1a'Vl 3 
Pp- 0 33 HP 



VauaUaa and conaUuLa l/om UaimokM wo<k*n*ai Unit 



r- 6#7 a 

c- 0 5411 ft 

b- 3 

A- 100 64 6*2 

dyn-aa 000133 

01* 1 407404 pa I 

BL - 0 033024 

Mo a 0 15003 



Wp- 18000 it* 

V4p- 204 02 Arc i a 

VU 40 00 knot* 

Afl- 10 44 

omaga 70 41 rad/a 

Cl a 0 00300 

Cl (avg)a 0 20 

la- 040 



*60 00 ip a 

47 11 Ip* 

780 i pm 



NACAO0I2 



CLmaa C4o Cl aipl-.a Mt< 

123 0 0047 0 73 0 71 



Aipn*l*JaM) 

13 



50 

59 

40 



70 

70 

40 

43 

00 

03 

100 



VH*P»> V«£it£*i 

07 03 480 00 

73 03 430 00 

44 30 450 00 

02 03 450 00 

101 27 480 00 

100 71 450 00 

110 13 450 00 

124 50 450 00 

133 03 450 00 

143 40 450 00 

131 00 450 00 

100 34 450 00 

104 70 450 00 



Pp IMP) Mu 

0 33 0 1500 

0 74 0 1400 

1 04 0 1075 

1 30 0 2043 

1 40 0 2230 

2 20 0 3430 

2 43 0 2023 

3 31 0 2013 

4 20 0 3001 

3 11 0 3100 

4 00 O3J70 

7 13 0 3543 

0 32 0 3751 



lam Cl 

0 0240 0 0031 

0 0243 0 0031 

0 0240 0 0031 

00203 0 0031 

0 0204 0 0031 

0 0313 0 0031 

00350 00031 

0 0104 0 0031 

0 0450 0 0031 

0 0013 0 0031 

0 0545 0 0031 

0 0440 0 0031 

0 0700 0 0031 



8 atyma 

0 07 JA 0 0013 

0 4730 0 0013 

00750 OOOI3 

0 0730 0 0413 

0 0730 0 0013 

0 0750 0 0413 

0 0750 0 0013 

0 0730 0 0013 

0 0730 0 0013 

0 0734 0 0013 

0 0734 0 0013 

0 0730 0 0013 

0 0730 0 0013 



0*(B '‘3*0 5*Mu ~ 2) 
33-IB-3* MUu'J'bl 
23* (B * 2*i8*' 3 ♦Mu'" 2W 
5*Uu*(8 *■ 2 ♦ 02a*Mu-2) 



All - (4*iMu*B * 2 13 3/S|)/(B ** 2(B - 3 0 5*Mu ~ 2)) 
AI2-(4*Mu-flJ/(3«(B“ 30 3*Uw~2}| 

A 13 a (2*Mu*B “ 2)/{B “ 3 0 5Mu *“ 3) 

All- IB *3* 1 8*Mvi ** 2)7(8'* 30 0*Mu “ 2) 



1 1-0 

12-0 

13- 0 

14- 0 

Li 12 

0 4/00 0 X<44 

0 4013 0 3003 

0 4020 0 3104 

0 4444 0 3114 

0 4440 0 3120 

0 4400 0 3143 

0 4414 0 3150 

0 4014 0 3174 

0 4007 0 3103 

0 4000 0 3311 

0 502* 0 3230 

0 5050 0 3231 

0 5003 0 3273 



T3 14 

0 2J01 0 0/ 10 

0 2314 0 0404 

02331 00407 

0 2340 0 0040 

0 2344 0 1041 

0 2340 0 1174 

0 2412 0 1207 

0 2434 0 1 342 

0 240 3 0 1404 

0 3400 0 1533 

0 2310 0 1040 

0 2340 0 1744 

0 2542 0 1044 



All A12 

1 1470 0 4154 

1 3321 0 4442 

1 2044 0 5233 

1 2013 0 3 770 

I 2243 0 0333 

I 1012 0 4402 

I 1544 0 7401 

• I 1314 0 4036 

I 0444 0 0424 

-1 0321 0 0223 

1 0173 0 0435 

0 0423 1 0457 

-0 0473 I 1004 



A 13 Al 4 

0 3037 I 0400 

0 3427 1 0010 

0 3422 I 0754 

0 4230 1 0010 

0 4424 I 1007 

0 5034 I 1204 

0 0440 I 1500 

0 5471 1 1743 

0 0300 1 1903 

0 0737 1 2303 

0 7143 1 2337 

0 7034 I 2070 

0 0102 1 3200 



l-o aquation* l**o unknown* 

3(Cl/*lgm*)/*al*mf 1 a lhai*0*12* thaiaTMl* I hai*2*l4 
tam*A1 1 ♦ Tha(*0*AI2* fhat*l»AI3* h>at*2*A 14 -0 
Solving kx T 1—1*0 and T ha La 2 

IliataO- ((haul *A 13* 1h*i*2*A 14 ♦ l*m*AI 1J/A 12 

1haU2- |A12*(2(Ci/*igm*)/a-t*m*T 1 lhaiat *T3> ♦ I2*l1-i*r *A 13 *lam*A1 1 *74/(14 *AI3 A14*t2} 



Ihaialaiwtti 

0 

• I £2t 

0 72 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 

0 73 0 



1hai*2 
-0 07 
0 07 
O 00 
006 
0 10 
-0 II 
0 13 
0 is 
0 10 
0 21 
0 25 
O 29 



T1 ataO *lpb*270 

0 00 0 00 

0 10 0 43 

0 10 0 00 

0 It 0 40 

0 13 >0 00 

0 13 13 23 

0 14 13 95 

0 10 10 11 

0 10 10 77 

0 22 22 04 

0 23 20 06 

0 30 30 07 



alpnaVO aipn*u70» - 1 1. a 1*0 fh*l*2* IhawT ♦ Lam/11 ♦ Mu) 
0 12 *1ph*{90) a JbaiaO* Ih*l*2 * lhalal *tam/(1 ♦ Mu) 
000 
000 
0 13 
0 20 
0 40 

0 73 
0 00 

1 31 
1 07 
200 
283 
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888888888888 8l£ 



Ihada RPV Maikop la. D*«| 



Pitch AngU D»tf mlnallon. 

Inflow >atio Um« (( TpW W/(iho*A*vn)/VMp) 

Advant* Ratio Mu • V(/V*p 

Paiadta Powa* (Pp)-O )Vho*W*VI * ) 

Pp- OU HP 



V*li*Na* and constant* bom Malnroto* link 



R- 8 #7 fl 

c- 0 54 17 « 

b- ) 

A- 10084 8*2 

dgma • 0 0111) 

0t- 1 *47 tM (,tl 

EU. • 0 031824 

Mu- 0 1900) 



Wg- 190 00 0>« 

VSp- 204 42 knoia 

VI- 40 00 knot* 

Afl- 10 44 

Orrwga 74 4 1 tad/* 

Cl* 0 00309 

CUavg). 0 20 

W- 080 



450 00 *>■ 

47 91 *>• 

759 i pm 



NACA0012 



V’ (*"<>'») 

40 4 7 )1 

45 75 99 

50 44 14 

98 9 2 4) 

80 101 27 

05 109 71 

70 119 18 

78 124 94 

40 1)9 02 

48 14) 48 

40 181 90 

49 140 54 

100 144 79 



CLmaa Cdo 

1 29 0 0097 

V»P <*P«) Pp (HP) 
450 00 0 9) 

450 00 0 74 

490 00 1 04 

450 00 IH 

450 00 1 40 

450 00 2 20 

450 00 2 8) 

450 00 ))1 

490 00 4 24 

490 00 8 11 

450 00 9 00 

450 00 7 1) 

490 00 4 )2 



CLatpba Mci 
5 7) 0 71 

Mu lam 

0 > 500 0 0244 

0 1844 0 024) 

0 1978 0 0249 

0 204) 0 0242 

0 2290 -0 0294 

0 24)4 00)1) 

0 2925 0 0390 

0 261) 0 0394 

0 )001 0 0450 

0)144 0 091) 

0 )379 0 0545 

0 354) 0 0440 

0 )781 00740 



Afpha{*l*ll) 

12 

Cl 8 

0 00)1 0 47)4 

00031 097)4 

00031 097)4 

000)1 097)4 

00031 017)4 

0 0031 0 9734 

0 00)1 0 97)4 

0 0031 0 97)4 

0 0031 0 97)4 

00031 0 97)4 

00031 0 97)4 

00031 0 97)6 

00031 097)4 



Thalal - Iwiat 



0 041) 
0 0913 
0 091) 
0 091) 



0 091) 
0 091) 
0091) 
0091) 
0 091) 
0 091) 
0 091) 
0 091) 
0 041) 



a I wial(iad) 

8 7) O 07 

5 7) -0 07 

9 73 -0 07 

8 7) -0 07 

) 7) 0 07 

9 7) -0 07 

8 73 O 07 

4 7) -0 07 

8 7) -007 

8 7) 4 07 

8 7) 0 07 

8 7) 0 07 

9 73 0 07 



1 1 «0 3*(B “2*0 5*Mu * 2) 
r].013*(B')i 05*Uu - 2*8) 
I)»0 29"(8 “ 2*(B * 2» Mu * 25) 

T 4 - 0 5 *Mu * (8 “ 2 ♦ 029 *Mu * 2) 



A 1 1 . |4‘|Mu*8 * 272 3,8)>/<B ~ 2(B “ 2 0 8*Uu * 2)| 
A 12- |8*Mu*8»7|)*(B* 2 0 )*Mu* 2)) 

A l). |2*Mu-8 * 2)/(B *24 5 Mu ~2) 

AI4. (B* 2*1 S*Mu * 2}/(B * 2 0 5*Mu * 2) 



T1 

0 4799 
0 441) 
0 4829 
0 4449 
0 4844 
0 4990 
0 4914 
0 4939 
0 4997 
0 4999 
0 9029 
09059 
0 509) 



12 
0 3084 
0)09) 
0 3104 
0 3119 
0 3129 
0)14) 
0 3184 
0 3174 
0)192 
0 3211 
0 3230 
0)231 
0 327) 



13 
0 2301 
0 2319 
0 23)1 
0 2349 
02384 
0 2389 
0 2412 
0 2419 
0 2482 
0 2489 
0 2914 
0 2949 
0 2582 



14 

0 0714 
0 0804 
0 0897 
0 0984 
0 1041 
0 1174 
0 1267 
0 1)42 
0 1454 
0 1852 
0 1944 
0 1744 
0 1944 



All 
1 3879 
1 3)21 
1 2984 
1 291) 
1 2292 
1 1912 
1 1984 
•I 1219 
1 0484 
-1 0921 
1 0172 
0 9823 
0 9473 



A12 
04194 
0 4892 
0 5232 
0 5779 
0 8332 
0 8892 
0 7491 
0 8039 
0 9829 
0 9228 

0 9835 

1 0457 
I 1094 



A 1) 
0 30)7 
0 34 27 
0 3822 
0 4220 
0 4824 
0 50)4 
0 9449 
0 5471 
0 9300 
0 97)7 
0 714) 
0 79)4 
0 9102 



AM 
1 0480 
I 0910 
1 0759 
1 0919 
1 1097 
1 1294 
1 1909 
1 1742 
I 199) 
1 2260 
1 2597 
I 2670 
I 3209 



Two aquation* two unknown* 

2(C Vdgrrw)/* • laml 1 » ThalaO* 12* Thai* T *13 ♦ 1hala2*T4 
Lam*A1 1 ♦ Tha(a0*A I 2 * Thalal *A13 * Th«U2*Al4-0 
Solving kx ThalaO and 7h«f«2 

1 h«i*o - (Thai.! *A I) ♦ Thau2*A 14 U am* A i |)/A 12 

Ihaia2- |A12*(2(Cl/dgma|/a Lam*f 1‘ThalaT*T3) *T2'ThaiaT*Al)a lam*A 11*1 2)7(1 4*AI2 AM*I2) 



Thai*) 
-0 07 
0 07 
009 
009 
O 10 
O 11 
0 1) 
0 15 
O 19 
0 21 
O 29 
0 30 



ThalaO *fpha270 
0 19 7 10 

0 15 7 47 

0 19 409 

0 14 484 

0 17 9 99 

0 14 II 40 

0 20 1) 17 

0 22 15 34 

0 24 19 II 

0 27 21 49 

0)1 29 99 

0 35 30 98 



a1ph*00 alpha (270) - ThalaO 1hai*2« ThalaT ♦Lam/(l ♦ Mu| 
0 92 alpha(90> • ThalaO ♦ 1haia2* IhataT ♦ l*m/(l ♦ Uu| 
O 95 
1 03 
1 1) 

I 30 
I 50 

1 7) 

2 01 

2 32 
209 
309 

3 99 
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|fi«t<* HPV rUucuv**' P-wgn 



Pilch Angl» 04 l»rmtn»llon. 

InNow i*mo l*in- i( f>V¥ W/(ino*A*V0)/V8p} 

UvviCI HiK) Uu-Vt/Wp 

Ptiuiu p»<mi (Fp)-o s*iho*u*vr 3 

Pp- OU HP 



V.H.U-* «nd cuniunu horn Mtvuakw wO'k*ri**l link 



H- 5 6/ 0 

c- 0 5417 H 

b- 3 

Am luo 64 n i 

0 00 133 

01- I 46/404 pU 

64.- 0 033624 

Uu « 0 15003 



My- 150 00 lb* 

V4p- 304 62 knot* 

VI- 40 00 knot* 

AA- 10 44 

om*y* 74 41 «*0/i 

Cl- 0 00300 

a (-»y|- 0 30 

t*- 0 50 



450 00 Ip* 
67 51 Ip* 

758 i pm 



NACA0012 




Clm»i Ctio Cl »lpti* Ucl 

1 25 0 008 7 5 7 3 0 71 



VI(’P») V»P 1»PH 

47 41 450 00 

73 05 *50 00 

44 34 450 00 

42 53 450 00 

■01 27 450 00 

108 71 400 00 

114 15 450 00 

124 30 430 00 

133 02 430 00 

143 44 430 00 

10 1 40 430 00 

100 34 450 00 

145 74 450 00 



Pptl'PI 

053 0 

0 74 0 

1 04 0 

1 34 0 

1 00 0 

2 26 0 

2 05 0 

3 51 0 

4 26 0 

Oil 0 

6 04 0 

7 13 0 

4 32 0 



1300 0 0244 

1444 0 0245 

1675 0 0244 

2043 O 0242 

2250 0 0244 

2434 0 0313 

2424 0 0350 

2D 1 3 0 0386 

3001 0 0450 

3164 4)0513 

3374 0 0645 

3543 0 0464 

3751 0 0740 



Cl B 

0 0031 0 8734 

0 0031 0 4734 

00031 06734 

0 0031 0 6734 

00031 06734 

0 0031 0 6734 

0 0031 0 6734 

0 0031 0 6734 

00031 06734 

0 0031 0 6734 

0 0031 0 6734 

00031 06734 

0 0031 0 4734 



In-u 1 . ttot*< 



00413 6 

00413 5 

00413 5 

0 0413 • 

00413 6 

0 0413 4 

0 0413 4 

0 0413 4 

00413 5 

0 0413 6 

0 0413 6 

0 0413 6 

0 0413 6 



• I »«l*l(l Xl) 

71 -0 14 

73 0 14 

73 6 14 

73 -0 14 

73 0 14 

73 O 14 

73 0 14 

73 A 14 

73 0 14 

73 0 1 4 

73 0 14 

73 A 14 

73 A 14 



11- 03*lB~2»03'Uu“2) 

12- 0 33*(B “3* 03 140 -*2-8} 
13.0 25’0B * 2*l8 '• 2* Mo ~ 2)) 

1 4 « 0 3*Mo*(8 ~ 2 1 025 14u “ 2) 



All- (4*|Mo*8 ** 2/2 ~ 2(8 * 20 5*Mo~2)) 

A 1 2- (6*Uv, *8 ) / (3 • (B ~ 2 0 5*Mo •* 2J ) 

A 13- (2*Mvi*B" 2|/(B ~ 2 0 3Uo * 2) 

A14- (8 ~ 2 i 1 3*Mu''2)/(B*20 314o * 2} 



M W 

0 4/86 0 X44 

0 4613 0 3083 

0 4624 0 3104 

0 «646 0 3116 

0 4344 0 3126 

0 4800 0 3143 

0 4614 0 3154 

0 4034 0 3174 

0 4847 0 3102 

0 4804 0 3211 

0 5024 0 3230 

0 5054 0 3251 

0 5043 0 3273 



12 11 

0 ^jOI 0 0/16 

0 2316 0 0404 

0 2331 0 0407 

0 2348 0 0868 

0 2344 0 1041 

02344 0 1174 

0 2412 0 1267 

0 2434 0 1342 

0 2462 0 M54 

0 2464 0 1552 

0 2516 0 1646 

0 2546 0 1746 

0 2562 0 1544 



All AI2 

1 *478 0 4154 

I 3321 0 4682 

1 2904 0 5 232 

1 2613 05770 

12202 0 6332 

I 1612 0 6482 

1 1544 0 7461 

t 1216 0 4038 

I 0446 0 6624 

10021 0 6225 

10172 0 6435 

0 0423 1 0457 

0 0473 11004 



A I J *14 

0 3037 I 0440 

0 3427 1 0410 

0 3422 1 0754 

0 4220 1 0014 

0 4624 1 1007 

0 5034 1 1204 

0 5446 I 1504 

05471 I 1742 
0 4300 I 1803 

0 4737 1 2245 

0 7103 1 2557 

0 7034 1 2070 

0 4102 13200 



1 wi »yu**o«i* two unknown* 

2(CV*lym*)/* -l*n\l I ♦ Ih*i*0’l2* lh»uf'U* It>*«*2*l4 
L*m*AI I ♦ lh«u0*A 12 * Ih*i*l *A13 ♦ Th*u2*A 14 -0 
Solving kM Th*l*0 *nd th*l*2 

11i*u0- (Ih*ul*Al3» Ili»u2* AI4 ♦ Imn’Al 1)/A 12 

l/>*<* 2- (A12‘{2|Cl/*40m*|/* LxnMl Th*l*I *T3) ♦ I2‘l1.*l*l *A| J »l*m*A1lM2)/|l4*A 12 A 14*1 2) 



lh»i*2 
0 07 
0 07 
004 
004 
0 10 
0 II 
0 13 
0 15 
0 IS 
0 21 
0 25 
0 30 



11—uO *lpb*2/0 

0 20 6 12 

0 20 4 51 

0 20 7 IS 

0 21 7 80 

0 22 9 12 

0 24 10 57 

0 25 12 39 

0 27 14 60 

0 30 17 45 

0 33 20 67 

0 37 25 07 

0 41 30 10 



*lp3i»90 *lph*(270| « 1h*u0 Ih*w2* lt.*Ul tUni/|l«Uu| 
■I 90 *lpA*(0O) • 1h*l*0* Ih*u2* 1 K*U f H*nV(HUu] 
1 06 
206 
2 17 
2 32 
2 52 

2 75 

3 02 
-3 33 

3 60 

4 10 
4 56 
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